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ABSTRACT 

A model describing the configuration of a water droplet evaporating on the 
surface of a radiantly heated semi-infinite solid is developed. A shape factor and the 
solid-liquid-vapor contact angle describe the transient droplet shape, though the initial 
value of the latter parameter is found to have a negligible effect on the droplet's 
evaporation. The droplet shape model and a modified radiation heat term are 
incorporated into a previously-developed computer model to predict the evaporation 
of a single droplet on a semi-infinite solid subjected to radiant heat input. The code 
predicts transient temperature profiles that agree well with experiment. A simplified, 
closed-form solution for the prediction of surface temperatures in the vicinity of an 
evaporating droplet is used to fit the data produced by the single droplet model. This 
closed-form solution facilitates calculations used in a model for the cooling of a 
surface by the evaporation of a sparse spray of water. The previously collected data 
base for sparse spray cooling using degassed water is expanded with a new set of 
experiments using water containing dissolved gases. Experimental results indicate 
that there is no significant difference in the sparse spray cooling when using either 
water containing dissolved gases or degassed water. Using superposition of effect, a 
computer code to model the sparse spray cooling of a solid surface subjected to 
radiant heat input from above is developed. The model agrees well with 
experimental data. The time constant governing the eciitine of the surface is found 
to depend on the properties of the solid surface and is independent of the water mass 
flux. For both experimental and theoretical cases, the solid used is macor. Macor is 


a ceramic with a low thermal-conductivity. The surface temperatures of interest 


iv 


range from 100 °C to the onset of the nucleate boiling condition for water. Radiant 
heat input is at relatively low temperatures (1.e., the radiant source is at less than 700 
i). 
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NOMENCLATURE 
radius of liquid droplet edge (see Figure 2) 
area 
specific heat 
droplet distribution function (see Equation 37) 
random droplet distribution factor (see Equation 39) 
frequency (see Equation 13) 
fractional coverage at polar angle @ (see Figure 8) 
interfacial radiative heat flux (see Equation 9) 


radiative heater coils to droplet surface view factor (see Equation 10) 


mass flux (see Equation 13) 

convective heat transfer coefficient 

Bessel functions 

thermal conductivity 

penetration depth (see Equation 15) 

heat flux 

average conductive heat flux (see Figure 14) 
total conductive heat energy term (see Equation 35) 
strength of point sink (see Equation 35) 

radial spatial coordinate 

radius of wetted region underneath the droplet 


droplet apex (see Figure 2) 
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time 
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temperature 
temperature drop 


droplet volume 
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Z axial spatial coordinate 
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a thermal diffusivity 
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6 thickness of macor solid 


€ absorptivity 

) solid-liquid-vapor contact angle (see Figure 1) 
© dimensionless temperature (see Equation 42) 
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total calculation area (see Equation 45) 
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1. INTRODUCTION AND BACKGROUND 

The removal of heat from a semi-infinite solid is the desired outcome for a 
number of engineering applications. Applying a sparse spray of water to its surface 
allows for large amounts of heat to be removed due to the high latent heat associated 
with the evaporation of water. Evaporative cooling and spray cooling are the 
common terms found to describe this kind of heat removal. 

Industrial uses for spray cooling include the quenching of molten metals 
during casting and the coating of surfaces to form protective finishes. Spray and 
mist cooling find a arity of uses in the power generation industry, such as the 
cooling of turbine blades and cooling tower applications. 

Several researchers have focused their studies on the fundamentals of the 
evaporation of droplets and their cooling effects. Simon and Hsu [1] studied the 
wetting characteristics of evaporating droplets on various surfaces. They recorded 
droplet shape histories at room temperature on copper, lucite and teflon surfaces. 
Both Toda [2] and Bonacina [3] performed early investigations of spray-surface 
interactions and provided fundamental insight into the uses of mist cooling. 
Photographic techniques were employed by Zhang [4] to determine flow patterns in 
evaporating droplets on glass and copper plates. Using these techniques, Zhang was 
able to construct liquid-air interfacial flow maps. 

In the area of fire suppression and protection, sparse spray cooling finds 
numerous uses. The current work represents the last phase in a research effort 


spanning several years to quanitfy and develop models for the spray cooling of hot 
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surfaces in a fire environment. The research has been developed and performed as a 
mutual effort between the Building and Fire Research Laboratory at the National 
Institute of Standards and Technology (NIST) in Gaithersburg, Maryland and the 
Mechanical Engineering Department of University of Maryland at College Park. 

Beginning in 1987, diMarzo and Evans [5] studied a single droplet 
evaporating on a high thermal conductivity surface. A eretiel model using a 
boundary element method to predict the cooling experienced by a semi-infinite solid 
due to an evaporating droplet was developed in 1989 by diMarzo, Kavoosi, and 
Klassen [6]. In 1992, Liao [7] developed a computer model to predict the transient 
thermal behavior of a solid caused by the evaporation of a single droplet and 
proposed a model for the impingement of a sparse spray of droplets. Both a high 
and a low thermal-conductivity surface heated from below by conduction were 
studied. Experimental techniques based on infrared thermography to record the 
evaporation of a droplet on a radiantly heated semi-infinite solid were developed in 
1992 by diMarzo and Kidder [8]. Dawson [9] extended this experimental work to 
record the effects of a random distribution of droplets (spray) on the surface. In both 
cases, the solid exhibited low thermal-conductivity. A stan model of the | 
evaporation of a single droplet for radiative heat input conditions constitutes the 
work of White [10] in 1993. ae 

With the progression as discussed above, the research described hereafter 
pertains to the development and validation of a computer model to predict the 


cooling of a semi-infinite solid subjected to a sparse spray of water where the heat 
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input is by radiation from above. Concurrently, a modification to the transient 
configuration of an evaporating droplet is proposed and an in-depth validation of the 
single droplet evaporation model for the case of radiant heat input against 
experimental data [9] is presented. Experimental work to determine the effects of 
dissolved gases in the water on the sparse spray cooling of a surface has also been 
performed and will be discussed. 


2. MODEL OF THE TRANSIENT DROPLET SHAPE 

The development of a computer model for the prediction of a single droplet 
evaporating on a hot solid surface requires consideration of the transient 
configuration of the droplet. The shape of the droplet can have a significant impact 
on its evaporation. A representation of the droplet geometry is determined for the 
case of a radiantly heated solid surface. This representation is incorporated into a 
computer model that predicts the evaporation of a single droplet on a semi-infinite 


solid subjected to radiant heat input [10]. 


2.1 Previous Geometric Model 
Common to all previous studies concerning the dropwise evaporative cooling 
of hot, solid surfaces [6,7] is the assumption that the droplet configuration is a 


spherical cap throughout its evaporation. This droplet geometry is given as 


(1 
a+y * (4+-y) 
Z= (— ) rt-— 


where, 


= = 
rR ty (oat +o (1425) a (2) 


8 is called the shape factor, or shape parameter, and is defined as the ratio of the 
radius of the wetted region under the deposited droplet to the radius of a spherical 
droplet of equivalent volume [3]. This early droplet model requires only a single 
parameter, namely the shape factor, to fully describe the droplet's shape as it 
evaporates. In this model, as the droplet evaporates, the radius of the wetted region 


under the droplet remains constant. 


2.2 Motivation for Developing New Droplet Shape Model 

A new model for the transient droplet shape is proposed in light of evidence 
found by both Chandra and Avedisian [11] and Zhang [4] which suggests that 
droplets take on a flattened configuration after deposition. The radiant heat field 
above the droplet is expected to affect the droplet shape in a manner simular to that 
suggested by Chandra and Avedisian. The increased temperature at the liquid-vapor 
interface caused by the radiant heat field above the droplet results in a relaxation of 
the droplet's surface tension. The relaxation causes the droplet to spread over a 
larger region of the solid surface and obtain a thinner, flatter configuration. 
Furthermore, the droplet model takes into account that upon reaching a given shape, 
the droplet's solid-liquid interface shrinks whereas it was previously assumed to 
remain constant throughout the droplet's evaporation. 

The shape of the droplet can have a significant impact on the evaporation 
process. Due to the thinner shape, a larger conductive heat flux from the solid to the 


liquid may occur. Therefore, it is important to accurately model the transient shape 
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of the droplet in order to capture its effects on the evaporation phenomena. 
2.3 General Description of the Droplet Shape Model 
As previously mentioned, the earlier droplet shape model requires a single 

parameter, namely the shape factor, to fully characterize its shape throughout 
evaporation. The model proposed here requires two parameters to fully describe its | 
shape at any time during evaporation. The first parameter is the shape factor, 8, 
while the second parameter is the solid-liquid-vapor contact angle, 0. 

The shape factor is determined from previous single droplet experiments [8]. 
The evaporation of the droplet basically occurs in two stages. The first stage is from 
the deposition of the droplet until the onset of the droplet's recession. Recession of 
the droplet is characterized by a shrinkage in the solid-liquid interface. The second 
stage is from the recession of the droplet until its complete evaporation. These two 
stages will be discussed in greater detail shortly. There are two values of the contact 
angle that need to be specified prior to using the droplet shape model. These are the 
initial value of the contact angle and the value of the contact angle upon the 
recession of the droplet. The initial value of the contact angle represents a Bein 
angle and the value of the angle at recession represents a minimum. The 


determination of these values are essential and will also be discussed. 


2.4 Determination of Initial Contact Angle Bounds 
The values of the initial and receding contact angles need to be specified 


prior to using the droplet model. The receding angle, 0, is the minimum solid- 
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liquid-vapor contact angle which results from a balance in surface tension and 
surface adhesion forces. The value of the receding angle was determined by Simon 
and Hsu using photographic experimental techniques [1] and confirmed by Kidder [8] 
to be 7°. The value for the initial contact angle is not as straightforward and in fact 
encompasses a range of possible values based on the initial shape of the droplet. In 
order to determine the range of initial contact angle values, the minimum and 
maximum limiting values of the range must be determined. 

The minimum value for the initial contact angle is found by assuming the 
droplet has a spherical Sine upon deposition. Therefore, knowing the value for 8, 


the minimum initial contact angle can be found by the following relationship 


aan 92. 1.(3) 


SimingatC val aay 


where ¥ is given by Equation (2). 

There are two possible values for the maximum initial contact angle 
corresponding to two situations which may occur upon deposition. First, a restriction 
is placed on the droplet shape such that the droplet has its maximum height upon 
deposition. This restriction requires that the initial height of the droplet be greater 
than or equal to the height of the droplet at recession. By setting the initial height of 
the droplet equal to the height of the droplet at the onset of recession, the 
corresponding maximum value for the initial contact angle can be determined. A 


second situation may arise, however, if the maximum initial contact angle 


corresponding to the recession height is greater than 90°. In this situation, the 
maximum value for the initial contact angle is taken as 90°, which agrees with values 
given by Chandra and Avedisian [11] for the solid surface temperature range of 
interest here. Figure 1 shows pictorially the two situations governing the maximum 


initial contact angle value. 


2.5 Transient Droplet Shape 

As the droplet evaporates, the transient droplet shape must transition smoothly 
from the initial configuration to its complete evaporation. The droplet shape 
transient consists of two stages, both of which provide smooth, gradual transitions 


between consecutive droplet configurations. 


2.5.1 Stage 1 - Initial Configuration to Onset of Recession | 

The first stage begins at the initial configuration and ends at the onset of 
recession. The initial configuration is determined uniquely by the initial contact 
angle, the droplet volume, and the radius of the wetted eon underneath the droplet. 
This radius is a function of the shape factor and is given by 


) 3 (4) 


If the contact angle is equal to its minimum limiting value, the initial 
configuration is that of a spherical cap shown at the top of the page in Figure 2 and 
described by the following equation 


In the model, s,' is iterated on until the resulting volume equals the initial volume. 


z=(p?=1?-b 


where, 


= R, 
tan (85, min) 


R 


aes fe} 
p sin (6 


Pentty, 


If the initial contact angle is not equal to its minimum limiting value, the 


(5) 


(6) 


configuration of the droplet can be described as a disk with a round edge as shown 


at the bottom of Figure 2. The round edge is assumed to fit a circular shape. The 


equations describing this configuration are 


i Seth 
Z=S'), TSR,-a 


mad Dee (Ale > Ras 


where, 


/ 
7a So 


oe (1+cos (6,) ) 


p 


/ 
3 So 


tT _ 
SCT (1+cos (8,) ) tan(—> -6,) 


(7) 


(8) 


The onset of recession is reached when the solid-liquid-vapor contact angle 
reaches the receding angle. At this point, the shape of the droplet is assumed to be a 
spherical cap and, replacing ©) ,,i, With 0, is described by (5) and (6). A spherical 
cap shape is expected because a minimization in the liquid-vapor interface for a 
given volume of water must occur before the droplet can recede. 

Once the initial configuration of the droplet and the configuration of the 
droplet at the onset of recession are known, a smooth transition between the two 
configurations can be determined. Referring to the top of Figure 3, two lines are 
drawn, one tangent to the droplet apex and one tangent to the solid-liquid-vapor 
contact point. The point of intersection of these two lines can be found for both the 
initial and receding configurations, labeled A and B in Figure 3. The line 
connecting points A and B represents the locus of intersection points for all 
intermediate droplet shapes. If the initial height of the droplet were equal to the 
height of the droplet at the onset of recession, the line connecting A and B would be 


horizontal. 


2.5.2 Stage 2 - Onset of Recession to Complete Evaporation 

Once the solid-liquid-vapor contact angle equals the receding angle, the 
droplet is assumed to have the shape of a spherical cap. Upon further volume 
reduction, the droplet radius and height decrease while the contact angle remains 


constant. As the droplet evaporates, the droplet aspect ratio remains constant. The 
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bottom drawing in Figure 3 illustrates this stage of the droplet's evaporation. 


2.6 Results Using the Proposed Droplet Shape Model 

The results of the droplet's evaporation using the new droplet model are 
shown in Figures 4 through 6. The figures presented here correspond to an initial 
surface temperature of 130 °C, an initial droplet volume of 9 microliters, and a shape 
factor equal to 2.3. Figure 4 contains curves of volume versus time using both the 
previous and the new droplet model, labelled A and B, respectively. For the new 
droplet model, various initial contact angles, ranging from a minimum value of 25° to 
a maximum value of 90°, are shown. The curves show that the droplet's evaporation 
time is lower for the previous model than for the new model. That is, the curve 
labelled A reaches a zero volume before those labelled B. Though the initial contact 
angle for the previous droplet model equals that of the minimum value for the new 
model, the recession of the solid-liquid interface that occurs in the new model affects 
the evaporation of the droplet. Therefore, it is found that the difference between the 
two models is ee after the droplet reaches its receding conditions. Figure 4 also 
shows that in the current model the value of the initial contact angle has very little 
effect on evaporation times. 

Other transient quantities that reflect the sensitivity of the evaporation process 
to the initial solid-liquid-vapor contact angle are the shape factor, B, the normalized 
radius, R/R,, and the solid-liquid-vapor contact angle, @. Figures 5 illustrates how 


these quantities behave for the minimum and maximum initial values of the contact 
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angle. Similar to evaporation time, the shape factor and the shrinkage of the wetted 
region (R/R,) are relatively independent of initial contact angle. The shrinkage of 
the wetted region (R/R,) does not display a smooth behavior due to the way in which 
the code redistributes the droplet volume after Sy time step, which will be 
discussed in more detail in Chapter 4. Another interesting result shown in Figure 5 
is that the receding angle is reached at the same time for the two limiting cases 
shown, at approximately 55 percent of the evaporation time. This agrees with data 
from diMarzo and Kidder [8] that suggests receding conditions are reached at 55 to 
85 percent of the evaporation time for initial surface temperatures between 130 and 
160,26. 

A final aspect worth investigating is the behavior of the transient heat flux 
under the droplet. In Figure 6, the heat flux is normalized with respect to the 
average heat flux transferred to the droplet through the solid-liquid interface during 
the entire evaporation process. In this figure, it is important to note that the largest 
flux, and thus vaporization, occurs at the droplet's edge. The shrinkage of the wetted 
region under the droplet is apparent at t/t = 0.9, where the radius of the wetted | 


region becomes approximately 40 percent of the initial value. 


2.7 Summary 

The modified transient droplet shape model presented consists of two 
parameters to fully describe its shape, as opposed to the single parameter 
characteristic of previous models. This new droplet shape model may take on a 


iz 


Chandra and Zhang [11, 4]. Results show that the evaporation time is longer for the 
droplet described by the new model than for the droplet described by previous 
models. The longer evaporation times were found to be a result of the droplet's 
recession. It was also shown that the new parameter describing the droplet 
configuration, namely the solid-liquid-vapor contact angle, only has an effect on the 
droplet's evaporation through its value upon. The transient droplet shape model 
discussed here is incorporated in a code [10] to predict the evaporation of a single 
droplet on a radiantly-heated, semi-infinite, solid surface. This code utilizes a 
boundary element method for the solid thermal behavior and three one-dimensional 
solution schemes for the liquid droplet. A predictor-corrector scheme couples the 
solid and liquid behavior. A validation of the computer model against experimental 


data is presented in Chapter 4. 
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Figure 1 - Droplet Configurations and Corresponding 
Maximum Contact Angles 
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Figure 2 - Initial Droplet Configuration 
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Figure 3 - Droplet Shape Transition 
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Figure 4 - Droplet Volume v. Time 
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Figure 5 - Effect of Contact Angle on Droplet Parameters 
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Figure 6 - Transient Normalized Interfacial Heat Flux 
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3. MODIFICATION OF THE RADIANT INTERFACIAL HEAT FLUX . 
IN THE SINGLE DROPLET EVAPORATION MODEL 


As mentioned earlier, a model was developed to predict the evaporation 
phenomena associated with a single droplet evaporating on a hot, low thermal 
conductivity surface where the heat input is by radiation from above [10]. A 
significant effort in developing this model centered around the quantification of the 
radiant heat input. While performing a recent validation of the code against 
experimental work, the author found that one of the parameters involved in the 
radiation portion of the model needed to be modified. The basis of that modification 


is discussed presently. 


3.1 Experimental Radiant Panels Geometry 

For the experimental work performed by Dawson [9] on the single droplet 
evaporation, heat input was by radiation from above the surface. The radiant heat 
input was provided by three, truncated conical radiant panels located above the 
surface. The geometry of the panels with respect to the surface is shown in Figure 
7. The two panels that are the same size have external diameters of 0.2 meters and 
depths of 0.8 meters. The third panel has an external pea of 0.3 meters and a 
depth of 0.07 meters and surrounds the perimeter of the surface. The panels can 
reach temperatures of 700 °C and are controlled by an Omega CN-7100 digital 


process controller. 
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3.2 Previous Determination of the Radiant Interfacial Heat Flux 

In developing the single droplet model to include radiation effects, White [10] 
divided the radiation into two parts with respect to the liquid, an interfacial heat flux 
term and a residual term that can be applied as a volumetric heat generation. White 


defines this interfacial heat flux as 
F,=F.oT,' (9) 


Viskanta and Toor [12] provide a general expression for F, which includes 
scattering within the liquid and bottom reflection. White has reduced this expression 


to 


F.=2f"£4cos (6) sin() (1-py) do (10) 


The terms making up this expression have the following meanings 


(1) - polar angle of incident radiation with respect to the normal to the 
droplet surface 


(2) - reflectivity from the droplet's surface, a function of 6 (less than 0.1 
‘for less than 60 °) 


(3) f, - fractional surface area coverage by the radiant panels (found by 
White for the laboratory geometry and shown in Figure 8) 


Several assumptions were made in determining a value for F,. These are 
(1) the liquid-vapor interface is horizontal and flat 
(2) scattering within the droplet is negligible 


(3) radiation reaching the liquid-solid interface is completely absorbed by the 
solid 


Sa 


(4) the geometry is symmetric, including azimuthal symmetry of the droplet 
and symmetry of the radiant panels with respect to the droplet 


(5) reflection at the liquid-vapor interface obeys Fresnel's relations for 
electro-magnetic theory - water is considered as a dielectric 


With these assumptions, Equation (10) was numerically integrated and a value of 
0.226 was calculated for F,. 


3.3 Modification to the Interfacial Radiant Heat Flux Term 

Of the assumptions listed at the end of section 3.2, the first one is the most 
questionable. Considering the droplet shape that is discussed in Chapter 2, the shape _ 
of the droplet surface, especially near the edges, may not have a normal that is in the 
vertical direction. If this is the case, much of the radiation coming in at large polar 
angles would not be reflected away as is assumed. 

To rectify the problem stemming from the curvature of the surface, ideally F, 
should be treated as a function of the droplet geometry. However, because 
evaluating F, in this way is cumbersome, a multiplier is instead introduced as an 
attempt to obtain an average value for F, that still includes the geometric effects. 

With the droplet shape considered as flat and horizontal, the term inside the 
integral given in (10) (f,cos(@)sin(@)(1-P,)) is evaluated at discrete values. The 
term is then evaluated for the droplet having a contact angle of 25° which is the 
minimum value discussed in Section 2.4. The angle, @, for this droplet geometry is 
actually the angle between the angle of incidence and the normal to the droplet 


surface, assuming the radiation passes through the center of the droplet. Taking a 


ee 


ratio of the quantity obtained after summation for the curved droplet geometry to the 
same quantity obtained for the flat droplet geometry yields a value of approximately 
1.9. Following the same procedure for a droplet having a contact angle equal to the 
receding angle yields a ratio of approximately 1. Therefore, the correction factor for 
F, should lie somewhere between 1 and 1.9. 

A modification of F, by a factor of 1.75 was used in the validation effort that 
is presented in Chapter 4. The average value between 1 and 1.9 of 1.45 was 
increased in order to provide a single, simple correction that incorporates other 


effects in the code that are not easily accounted for. These other effects will be 


expanded upon in Chapter 4. 


oe 
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Figure 7 - Experimental Radiant Panels Geometry 


24 


| o LL 
| pal 
| | © 
: | | es 
| I = 
Sain 
O 
| ab 
pefenboped 
= ne oY if aa 
S s S 
SJOVYSAOO WNOILOVYS 


‘Figure 8 - Fractional Surface Area Coverage, f, 
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4. VALIDATION OF THE SINGLE DROPLET COMPUTER MODEL 

In a previous model of the sparse spray cooling of semi-infinite solids [7], a 
database was used to provide the temperature distribution during the vaporization of 
a droplet. The paramters of the database consisted of various times during the 
evaporation of a single droplet, various initial surface temperatures, and the 
corresponding temperature drop. The information comprising the database was 
obtained from a single droplet computer model. Both the single droplet model and 
sparse spray model in the previous work focused on a solid that was heated by 
conduction from below. 

The sparse spray cooling model of current concern is discussed in detail in 
the Chapter 6 and focuses on the cooling of a surface that is being heated from 
above by radiation. Again, in order to accurately model what is happening at a point 
on the surface relatively close to an evaporating droplet, the single droplet model 
mentioned in the two preceding chapters needed to be studied and quantified. 
Results of the validation of the single droplet model are now discussed. A closed- 
form solution that approximates the single droplet's effect on the surface will be 


presented in section 4.5. 


4.1 Brief Description of the Single Droplet Computer Model 

The computer code of White [10] models a single droplet evaporating on a 
surface that is under the experimental conditions shown in Figure 9. The results 
used as the basis for validation of the sparse spray cooling model were also obtained 


using this experimental setup. This section is intended to provide a general idea of 
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how the single droplet computer code works. Details of the development of the 
model can be found in reference [10]. 

| The droplet shape is modelled using the geometry proposed in Chapter 2. 
Radiation effects are handled in the manner described in Chapter 3, with the largest 
portion of the radiation applied as an interfacial heat flux and the remaining 
distributed throughout the droplet as a heat generation term. The modification to the 
interfacial heat flux term discussed in Chapter 3 is used in the validation presented 
here. 

The model assumes a one-dimensional conductive heat transfer with internal 
generation in the liquid layer. The one-dimensional assumption was justified 
previously from diMarzo's coupled model [13] and is reinforced by the thinner 
droplet seomnetty that results from the direct radiation. The thermal behavior of the 
solid is governed by a transient, two-dimensional (r,z) diffusion equation which is 
solved using a Boundary Element Method. Coupling of the liquid and the solid is 


accomplished using a predict-correct method. 


4.2 Validation of the Single Droplet Computer Model 

Experimental results obtained by Dawson [9] for the evaporation of a single 
droplet on a macor surface are used as the basis for validation of the single droplet 
computer model. Validation of the code is determined by two criteria. The first 
criterion is the total evaporation time of the droplet. The second criterion is the 


transient surface temperature profile predicted by the code. 
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4.2.1 Inputs to the Computer Model 

Running the computer model requires several inputs. These inputs are used 
as a description of the conditions that occurred during an experiment. The first two 
parameters that are discussed relate to the initial heating conditions of the 
experiment, that is the initial surface temperature and the temperature of the radiant 
panels. The single droplet evaporation experiments run by Dawson were for initial 
surface temperatures of 110 °C and 130 °C. Therefore, these are the initial surface 
temperatures that are used for validation purposes. Various initial surface 
temperatures correspond to radiant panel temperatures when a temperature of the 
solid 0.0254 meters below the surface is held constant. The experimental surface 
bottom temperature ranged from 32 °C to 40°C for the various initial surface 
temperatures and was maintained by a chill plate underneath the surface. A curve fit 
to the panel temperature versus initial surface temperature was obtained from 
experimental data and is shown in Figure 10. Using the curve fit, panel temperatures 
of 473 °C and 512 °C are used for initial surface temperatures of 110 °C and 130 °C, 
respectively. 

Three remaining parameters are needed in order to use the model. These 
parameters pertain to the droplet itself and are the initial droplet volume, the initial 
shape factor, and the initial and receding contact angles. Dawson measured an initial 
droplet volume of 10 microliters for his single droplet experiments. For the 
validation of the code, an initial droplet volume of 9 microliters is used. This 


volume is selected to represent the average droplet volume that was used during the 
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experiments. A volume decrease is expected to occur because an experimental run 
consists of consecutive droplet evaporations and surface recoveries. In this way, 
while the surface recovers to its initial temperature, the water in the dispenser is 
being heated. Some of the water will evaporate, thereby causing less water to be 
dispensed for the following cycle. 

The initial shape factor, B,, is the ratio of the diameter of the droplet 
immediately after impact to the diameter of the equivalent spherical droplet prior to 
impact. For initial surface temperatures ranging from 110 °C to 160 °C, Kidder 
found initial shape factor values in a range between 1.5 and 2.5 [8]. The initial 
surface temperatures of interest here are 110 °C and 130 °C, yielding shape factor 
values approximately equal to 2. A value of 2.3 is used for the code validation as an 
attempt to account for the difference in droplet release height between the 
experiments of Kidder and Dawson. Dawson's release height was 0.3 meters higher 
than Kidder’s. 

The geometric model presented in Chapter 2 requires a value for the initial 
contact angle in addition to a value for B,. Because the study of that model showed 
that the initial contact angle had no effect on the droplet's evaporation, a minimum 
value for the contact angle is selected for the validation purposes. The minimum 
contact angle corresponding to the values for B, and the initial droplet volume is 
found to be approximately 25 °. Therefore, the droplet shape model that is employed 


is that of a spherical cap with the recession of the solid-liquid interface. 
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4.2.2 Validation Results 

Using the various input parameters stated in 4.2.1, evaporation times and 
surface temperature distributions were obtained. Evaporation times predicted by the 
computer model were 45 seconds for an initial surface temperature of 110 °C and 26 
seconds for an initial surface temperature of 130 °C. The experiments selected for 
comparison had evaporation times of 40.2 seconds and 24 seconds, respectively. 

The difference in evaporation times results from the discretization of both the. 
time and the water droplet in the computer model. The water droplet is modeled as 
discrete annular columns in the code. During each time step, the amount of water 
that makes up a particular column is the maximum amount of water that can be 
vaporized in that column, even if the heat supplied to that column during a time step 
is more than is required to vaporize the column. For the next time step, the 
remaining droplet volume is redistributed and the columns are reconfigured. In 
reality, the time and droplet volume are continuous so that all of the heat that is 
supplied acts to vaporize the droplet. Part of the redistribution effect is corrected 
through the increase in the value of F,. An error in the evaporation time on the 
order of 10% is acceptable in view of the arbitrary correction. Generally, errors are 
smaller than 10% except at low temperatures where the vaporization process is 
longer and redistribution effects are larger. 

The predicted surface temperature distributions are compared to the 
Seoeal distributions in Figures 11 and 12. The raw data used in these figures 
are contained in Table 1 and Table 2. The figures are for initial surface temperatures 
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of 110 °C and 130 °C. Three times during the droplet's evaporation are displayed. 
The first time corresponds to the beginning of the droplet's evaporation, the second 
time is just prior to the complete evaporation of the droplet, and the last time is just 
after evaporation of the droplet when the surface is recovering. 

A single line depicts the predicted temperature distribution, while the thicker, 
jagged line represents the experimental data. The appearance of the experimental 
results are due to the means of data acquisition. A line scan function on an infrared 
camera is used to display the temperature distribution. Using a PC and frame 
grabber software, the line is digitized in a pixel by pixel manner. Each pixel 
location in the vertical direction can be associated with a temperature and each pixel 
in the horizontal direction can be associated with a radial location. The thickness of 
the data arises because the line scan function produces a band of pixels Bien it is 
grabbed and digitized. Furthermore, note that during the droplet's evaporation, 
temperature distributions begin at the edge of the droplet. Using the thermographic 
technique, recorded surface temperatures underneath the droplet cannot be associated 
to infrared emmissions. 

From the comparisons made in Figures 11 and 12, it is evident that the 
computer model is able to predict reasonable transient surface temperature 
distributions caused by the evaporation of a single droplet. At early times, there is 
some deviation between the predicted and experimental temperatures at the droplet's 
edge. Again, this is a result of the inability of the need camera to record what is 
happening underneath the droplet and to the experimental uncertainties associated 
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with the position of the droplet's outer edge. As the surface recovers, the computer 
model lags the experimental data slightly for normalized radial location between 0 


and 1. However, as time progresses, these differences are decreased. 


4.3 Closed-Form Solution 

A closed-form solution for the transient surface temperature distribution is 
now presented. The closed-form solution provides a simple, yet reasonable means of 
obtaining temperature differences at points on a surface near the droplet. It is useful 
for the model of sparse spray cooling because it eliminates the need for a database, 


which can take up significant amounts of computer memory and computation time. 


4.3.1 Description of the Closed-Form Solution 
A closed-form solution to the problem of a disk with constant and uniform 
heat flux on an insulated surface is given in reference [14]. The solution has the 


following form 


Tease fala ae erf | SE gh qa) 
where q is the constant and uniform heat flux acting over a disk of radius R. The 
temperature at various locations (r) and times (t) can be found using Equation (11). 
Although the evaporation of a single droplet does not have a constant and 
uniform solid-liquid interfacial heat flux, the solution described by (11) models the 
phenomenon closely if minor modifications are made. Specifically, a modification to 


the radius, R, needs to be made in order to account for both the recession of the 
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solid-liquid interface and the fact that the conductive heat flux is neither constant 
nor uniform. The heat flux, q, is replaced by the term (q, + q,), where q, accounts 
for the conductive contribution to the vaporization of the droplet and q, is the initial 
steady state heat flux through the solid. Equation (11) yields a solution for the 
surface temperature distribution during the evaporation of the droplet. An expression 
can also be found for the solid surface temperature distribution for times after the 
droplet has evaporated (t > T) by modifying (11) to obtain [10] 


eat ee ane ea) teri tt) ere Veet) 1 a (12) 
Rilo altar R R x 


Surface temperature distributions in the vicinity of a droplet can be obtained for the 
evaporation of the droplet as well as for the recovery of the surface using Equations 
(11) and (12). However, the conductive heat flux contribution to the droplet's 
evaporation (q,) and the evaporation time (T) must be known before the solutions can 


be used. 


4.3.2 Determination of q, and T 

In order to use the solutions given in Equations (11) and (12), the evaporation 
time (T) and the conductive heat flux contribution (q,) to the vaporization of the 
droplet must be known. The initial heat flux (q,) is known from the initial surface 
temperature, the bottom surface temperature, and the thickness of the solid. The 
evaporation time can be obtained experimentally or by using the single droplet 
computer model. The conductive heat flux must be found using the computer model. 


For a range of initial surface temperatures and their corresponding panel 


SS 


temperatures, with the other inputs discussed in section 4.2.1, the computer model 
was run to obtain the resulting evaporation times and conductive heat fluxes. Results 
are plotted in Figure 13 and Figure 14, and curves are fit to the computed data 
points. The results are used in the comparison of the closed-form solution against 
the computer model. 


4.3.3 Temperature Distribution Results 

Surface temperature distributions given by the closed-form solution are 
compared against those obtained using the computer model in Figures 15 and 16 and 
raw data are in Tables 3 an 4. The distributions shown are for the same conditions 
used for the validation results of the single droplet model. Comparisons are made 
for two times during the droplet's evaporation and one time after the droplet has 
evaporated for initial surface temperatures of 150 °C and 100 °C. These initial 
surface temperatures represent a relatively high and a relatively low value for the 
range of surface temperatures that are found in the spray cooling database of concern 
to fire protection applications. 

In order to compare the closed-form solution to the results of the single _ 
droplet model, the radial shrinkage of the droplet needs to be accounted for. Since 
the solution represents a fit to the data, a simple means to account for the shrinkage 
of the wetted region is to introduce a constant into the expression. The constant 
effectively decreases the temperature difference given by the closed-form solution. 
Such a decrease has the same effect of shrinking the radius. The effect of 
introducing this factor into the closed-form solution is most apparent in the first two 


34 


temperature distributions shown in Figure 15. The first distribution predicted by the 
closed-form solution lies slightly above that predicted by the code because of the 
modification to the radius. The next distribution, which corresponds to a time after 
the recession of the droplet has begun, predicted by the closed-form solution lies 
slightly below that predicted by the code. If the modification to the radius was not 
included, both of the closed-form distributions would lie below the code distributions 
with the latter one significantly lower. 

Temperature distributions predicted using the single droplet computer model 
are shown by thick, solid lines and those obtained from the closed-form solution are 
shown by the thin, dashed lines. The closed-form solution for the temperature 
distribution agrees well with the predicted distributions with respect to both space 
and time. In these figures, distributions underneath the droplet are shown. The 
difference in character between the closed-form solution and the single droplet model 
underneath the droplet is a reflection of the boundary conditions that govern the two 
problems. The closed-form solution has a constant heat flux boundary condition 
underneath the droplet, while the computer model has a constant temperature 
undemeath the droplet. The differences in temperature distribution underneath a 
droplet are not expected to cause significant error in the sparse spray cooling model. 
Since a radial temperature distribution corresponds to the droplet, the accuracy of 
predictions near the droplet is less important than the accuracy far from the droplet 
due to the area associated with each. Additionally, experimental data does not yield 
correct temperature distributions underneath a droplet. 
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At an early time during the droplet's evaporation, there is virtually no 
difference in the distributions except for the behavior underneath the droplet which 
was just discussed. The biggest difference between the distributions is found at a 
time just prior to the complete evaporation of the droplet. At this time, the shrinkage 
of the droplet has occured and is apparent in the distribution predicted by the single 
droplet model. The effect of using the closed-form solution as opposed to the single 
droplet model to calculate the average surface temperature in the sparse spray 
cooling code will be quantified in Chapter 6. Immediately after the droplet has 
evaporated, again slight discrepancies are shown between the two curves undemeath 
the droplet. However, these discrepancies die out quickly as time progresses. 


4.4 Summary 

The single droplet computer model of White was validated with respect to 
evaporation times and surface aaneraine distributuions against experimental data. 
A relatively simple closed-form solution was shown to yield temperature distributions 
that agree with the model predictions. The closed-form solution requires the droplet — 
evaporation time and the conductive heat flux contribution to the vaporization of the 
droplet. These two pieces of information can be obtained using the single droplet 
model. A model for the sparse spray cooling of a surface for the same experimental 
setup as described in this chapter utilizes the closed-form solution as a simple means 
of calculating temperature depressions at points in the vicinity of droplets. The 


sparse spray cooling model is the subject of Chapter 6. 
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Figure 9 - Experimental Setup 
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Figure 10 - Panel Temperature v. Initial Surface Temperature 
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Figure 11 - Surface Temperature Distributions; 
T,, = 110°C 
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Figure 12 - Surface Temperature Distributions: 
T,, = 130 °C 
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Table 1 - Raw Temperature Data for Figure 11 
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Table 2 - Raw Temperature Data for Figure 12 


42 


100 


eo 80 
aw = 1300*exp(-0.03*Tso) 
iD) 
E 60 
- 
2 
m 40 
; 
© 
oe 
> 20 
jaa 
0 : 
90 100 110 120 130 140 150 


Initial Surface Temperature (C) 
Figure 13 - Evaporation Time v. Initial Surface Temperature 


—_-_ NY NY WO WW 
Aa OS th ae 
Go oo &G G&G & 
See =|] |} = 
oO ©& Oo) oF 


10000 


Average Conductive Flux (W/m?) 


17) 
i=) 
i=) 
cn) 


90 100 110 120 130 140 150 
Initial Surface Temperature (C) 


Figure 14 - Average Conductive Flux v. Initial Surface Temperature 
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Figure 16 - Surface Temperature Distributions; 
100 °C 
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Table 3 - Raw Temperature Data for Figure 15 
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Table 4 - Raw Temperature Data for Figure 16 
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5. EFFECTS OF DISSOLVED GASES ON THE SPARSE SPRAY 
COOLING OF SURFACES 


An experimental study performed previously examined the evaporative 
cooling of a radiantly heated, low thermal conductivity material subjected to a sparse 
spray of deionized and degassed water [9]. The transient behavior of the surface 
with respect to the average temperature and the spatial surface temperature 
distribution was quantified and is used in Chapter 6 to validate the sparse spray 
model results. The experimental work is extended here to include the effects of 
dissolved gases in the water on the transient thermal behavior of the surface. 


5.1 Description of the Experiment 
The experimental setup and procedures basically follow those described in 
more detail by Dawson and diMarzo [9]. A brief description of the experimental 


setup and procedures is provided here. 


5.1.1 Experimental Setup 

The schematic of the experimental setup is shown in Figure 9 of Chapter 4. 
The setup consists of a ceramic macor tile that is heated by three conical radiant 
panels from above and maintained at a constant temperature below by means of a 
chill plate. Two of the radiant panels are the same size and mounted above the 
macor surface at angles of 30° from the vertical on either side. The third panel is 
larger in diameter and shallower in depth and surrounds the perimeter of the tile. 


The temperature of the panels are controlled by an Omega CN-7100 digital process 
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controller and can achieve a maximum temperature of 800 °C. 

A droplet dispenser, shown in detail in Figure 17, hangs above the surface 
and, together with a set of solenoidal-controlled bumpers, provides a random spray of 
water to the surface. The dispenser was designed and tested previously [15]. The 
dispenser consists of a tapered, conical, aluminum body with a bored out center. The 
cavity exits through a hole at the bottom of the main body. This cavity is filled with 
water and sealed at the top with a plastic diaphragm and o-ring. A size 20 IV needle 
screws into the hole at the bottom of the cavity. A steel piston rests on top of the 
diaphragm while a SOIR acet mechanism is fitted to the top of the piston. The 
spacer provides a gap, which can be adjusted by a small screw, between the solenoid 
core and solenoid body. When the solenoid is energized, it pushes down on the 
spacer, closing the gap and causing the piston to deflect the diaphragm and thus 
produce a droplet from the needle. The size of the droplet is controlled by 
adjustment of the gap and the frequency of droplet deposition is controlled by the 
solenoid. 

The cavity in the dispenser is kept filled by a feedwater line which, at one 
end, fits into a small inlet located just below the diaphragm. The other end 
connects to a reservoir of water that produces a static head of almost zero. The 
height of this reservoir is adjusted throughout the experiment to maintain the near 
zero static head. 

An infrared camera looks through a chill pipe that ends just above the macor 


tile. The chill pipe is painted black on the inside and serves as a means of absorbing 
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any stray reflections from surrounding objects. The infrared camera can then focus 


on the surface of interest and record its thermal variations. 


5.1.2 Experimental Procedure 


A series of experiments were run for various initial surface temperatures 


ranging from approximately 110 °C to 180 °C. At each initial surface temperature in 


an experimental session, three mass fluxes of the gassed water impinging on the 


surface were tested. The initial surface temperature and mass fluxes were selected to 


closely match those from Dawson's experiments with degassed water. The procedure 


that was followed for each experimental session is outlined below. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


The macor surface is cleaned with ethyl alcohol and a soft cloth, 
lightly rinsed with distilled water, and allowed to dry. 


The radiant heaters are turned on and allowed to heat up for 
approximately two hours. 


The infrared camera is turned on, also two hours prior to 
experimentation. 


The reservoir is filled with distilled water and adjusted to create 
a zero static head so that no water flows from the IV needle. 


The droplet generator is turned on and allowed 10 - 15 minutes 
to stabilize. Once stabilized, 50 droplets are collected in a 
beaker which is quickly capped to avoid evaporation. The 
beaker and drops are weighed and the volume of a single droplet 
is determined. | 

Once the droplet size is adjusted and determined, the frequency 
of the solenoid motion is set to the desired mass flux. The mass 
flux is given by 


Gab VE (13) 
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where p is the density of water, V is the droplet volume, f is 
the frequency of droplet impingement, and A,, is the droplet 
impingement area determined by Dawson to be 0.0033 m’. 

(7) Once the surface temperature has reached the desired, steady 
value (as indicated by an Omega thermocouple probe), the 
droplet dispenser is set in motion with the solenoidal bumpers 
and allowed to stabilize. 

(8) | After the above procedures are completed, experiments begin. 

. Droplets are allowed to fall for about twenty-five minutes. Each 
run is recorded by the infrared camera onto 8 mm videotapes 
from which the data is later processed. 

5.1.3 Data Acquisition and Reduction 

The transient thermal behavior of the surface is captured using infrared 
thermographic techniques that were developed and implemented previously by Klassen 
[6]. An Inframetrics Model 525 infrared camera detects radiation wavelengths from 8 
microns to 12 microns and translates thermal variations of an object into a real time 
gray image. These images are made up of dark shades that represent cool regions and 
light shades that represent hot regions. The camera records the thermal image of the 
surface onto 8 mm videotape for future processing and convenient storage. 

The means of reducing the stored data was developed and used by Dawson [9] 
for his multi-droplet evaporative cooling experiments. A Matrox MVP-AT frame 
grabber board is installed in an IBM PC-AT and used to digitize single frames of the 
filmed images. Once digitized, each frame can be analyzed pixel by pixel using 
Imager-AT software. Each frame is calibrated using a temperature versus intensity 
relationship so that shades of gray may be translated into corresponding temperatures. 
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For the purposes of the multi-droplet evaporative cooling experiments, the 
infrared camera records a region of the macor surface of approximately 0.075 m x 
0.055 m. The pixels contained in this viewing area are 512 by 480 in the horizontal 
and vertical direction, respectively. Thus, the average spatial resolution is 0.00013 
meters per pixel. There are 130 shades of gray to represent an average temperature 
range of 100 °C. The temperature resolution is, therefore, calculated as 0.77 °C / gray 
shade. 


5.2 Results 

The results of both the dissolved gases and degassed experiments are presented 
and compared in Figures 19 through 32 and in Tables 6 through 12. The riansiens 
thermal behavior of the surface is represented in the form of average surface 


temperature versus time plots and in the form of surface temperature contour plots. 


5.2.1 Transient Average Surface Temperature 

Figures 19 through 23 contain plots of transient average surface temperatures. 
These plots were made using the computer program, PIXAV [9]. PLXAV is used to 
average across 3717 pixel values contained in an area of 0.046 m x 0.034 m and 
provide the corresponding average temperature. The average surface temperature is 
obtained in this way at 30-second intervals throughout an experimental recording. 

Seven experiments are compared. These experiments consist of initial surface 
temperatures of 111 °C , 131 °C, 151 °C, 162 °C, and 182 °C. The lowest three 
temperatures correspond to evaporative cooling, while the highest two temperatures 
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correspond to the onset of nucleate boiling and full nucleate boiling, respectively. At 
least one mass flux is compared at each initial surface temperature with an additional 
mass flux for initial temperatures of 151 °C and 162 °C. 

Each plot contains both raw data points from the dissolved gases and degassed 
experiments and a best fit line corresponding to each set of raw data. The general 
trend apparent in each plot is the decay of the average surface temperature from its 
initial value to some steady state value. The deviation of the data points from a 
smooth decay occurs due to the nature of the data acquisition. Because only a portion 
of the sprayed area is viewed and averaged, at any instant the number of droplets that 
can be seen may be different than at other instants thus resulting in fluctuations of the 
average surface temperature. 

The method used to determine the best fit curve to the data is the same method 
used by Dawson [9]. An equation of the following form is used as a curve fit to the 
temperature data | 

T= (T)-Ty,) © +7 yg (14) 
This equation is basically a decaying exponential with T, = initial surface temperature, 
T,, = steady state surface temperature, and T = time constant. The steady state 
temperature and time constant are adjusted so that an equal number of data points fall 
above and below the curve. 

Dawson suggests that the time constant is a function of the material as opposed 


to a function of mass flux [9]. Table 5 lists the time constants corresponding to the 
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various initial surface temperatures and mass fluxes for both sets of experiments. The 
time constants listed are fairly constant ranging from 126 to 270 seconds, where the 
average value is 186 seconds with a standard deviation of 11.5 seconds. The 
unchanging time constant supports Dawson's conclusion of the material dependence of 


this parameter. 


Table 5 - Time Constants for Sparse Spray Cooling Experiments 
(* indicates dissolved gases experiments) 


Using the limiting values found for the time constant, the penetration depth, 1, of the 


solid can be found from 
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T (15) 


where, 0, is the thermal diffusivity and has a value of 5.79 * 10” m/s for macor. The 
penetration depth is then calculated to be between 0.85 and 1.25 cm. The radius of 
influence of the droplet on the solid surface is a circular region with a radius of 3.5 to 
5 times the deposited droplet's radius [16]. For a 9 microliter droplet with a shape 
factor of 2.3, which are the conditions used for the single droplet validation, the radius 
of influence is approximately 1.05 cm to 1.5 cm. These values are of the same order 
of magnitude as the range above and, again, indicate that the time constant is a 
function of the solid material. 

Overall, the plots indicate that dissolved gases in the water do not affect the 
cooling of the surface. This is indicated by the temperature correlations to the best fit 
curves. For comparable initial surface temperatures and mass fluxes, the best fit 
equations to the data are similar in all cases. In light of the relative independence of 
the time constant discussed previously, the steady state temperature becomes the 
parameter of interest for comparison purposes. The average difference in steady state 
temperatures for the cases presented is 3.4 °C. The largest difference is 8 °C and is 
found at initial surface temperatures of 111 °C and 131 °C. Part of this difference may 
be explained by the difference in mass flux used for the dissolved gases experiment, 
thereby altering the steady state temperature. Additionally, the droplet distribution 
may not be the same in both experiments and therefore calculating an average surface 
temperature over only a portion of the entire surface could bias the results. 
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5.2.2 Surface Temperature Contours 

Spatial surface temperature results of the dissolved gases and degassed 
experiments can be found in Figures 24 through 32. These results are in the form of 
surface temperature contour plots. 

Three of the cases that are contained in Figures 24 through 30 are shown in 
order to convey the overall qualitative behavior of the surface. The cases shown are 
for initial surface temperatures of 111°C, 151°C, and 162 °C with mass fluxes of 
approximately 0.2 g/m’s, 0.96 g/m’s, and 0.97 g/m’s, respectively. 

Each figure contains two contour plots that correspond to the dissolved gases 
and degassed experiments. Three different times during an experimental run are 
selected. The times are chosen to show 1) a very early time in the transient when the 
surface temperature is close to its initial value, 2) a time fairly far into the transient, 
but before steady state has been reached, and 3) a time after steady state has been 
reached. , 

The contour plots for both the dissolved gases and degassed experimental 
results contain some features that are similar. Both sets of results show very distinct . 
locations where droplets are evaporating on the surface, or have just evaporated from 
the surface. These places are marked by concentric contours that increase in 
temperature as the outer contour is approached. Places where droplets have recently 
evaporated will still have concentric contours but the inner contour will be associated 
with a temperature higher than that of evaporating droplets. At a given time, roughly 


the same number of droplet locations can be seen on each contour plot. There may be 
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slightly more or less drops shown due to the fact that some of the droplets fall outside 
of the area showed by the contour plots. Another similarity between the contour plots 
for experiments is the way in which the plots progress with time. At later times, the 
effects on the surface temperature due to individual droplets tend to merge. While at 
earlier times, temperature effects are contained within a local region around the 
droplet. Lower temperature contours are found near the perimeter of ) the plot at later 
times than those found at early times, indicating that the entire surface is cooling. By 
a time of 600 seconds, the average surface temperature, which is labelled at the top of 
each plot, has approached the steady state temperature. | 

Alternatively, there are features of the contour plots that distinguish the 
behavior of the water containing dissolved gases from the behavior of the degassed 
water. The most noticeable difference between the two sets of results is the size of 
the concentric contours surrounding a droplet. Experimental results using water with 
dissolved gases show contours that are larger in the region surrounding a droplet than 
those corresponding to degassed water experimental results. Furthermore, the 
temperature gradients surrounding the water droplets containing dissolved gases are 
also larger. Large temperature pen are illustrated by the proximity of two 
concentric contours. 

The explanation for these differences stems from the make up of the two types 
of water droplets. Figure 18 depicts a general concept of the two types of water 
droplets and the behavior of the conductive heat flux associated with them. Note that 
no study was performed on the actual make up and behavior of a single droplet, 
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rather Figure 18 represents what the author believes to be occurring. The larger 
contours associated with a droplet of water containing dissolved gases can be 
attributed to the spreading of the droplet on the surface due to the gases. For droplets 
of equivalent mass, a droplet containing dissolved gases requires more area on the 
surface due to the relatively larger volume taken up by the gas. Consequently, the 
dissolved gases push the liquid outward. Higher temperature gradients are also 
explained by the presence of the gases. In locations where gas bubbles exist, a 
degradation in heat flux through the droplet occurs. Significantly larger heat fluxes 
are therefore found at the droplet's edge. These large heat fluxes cause high 
temperature gradients at the edge of the droplet. In contrast, the heat flux through a 
degassed water droplet is fairly constant with a smaller increase at the edge. Most 
likely, overall surface temperatures do not differ significantly for the two experiments 


because the average heat fluxes associated with the droplets are similar. 


5.4 Summary 

Transient average surface temperatures and surface temperature contours were 
compared for water containing dissolved gases and degassed water. The results 
showed that there is no significant effect on the cooling of the surface when using 
water containing dissolved gases as opposed to degassed water. The characteristics of 
surface temperature contours differ for the two types of water. The spreading of the 
droplet on the surface and the high temperature gradients at the droplet's edge result 


from the dissolved gases in the water and are portrayed in the contour plots. 


58 


The work presented in this chapter does not represent an in depth study of the 
evaporative behavior of water containing dissolved gases, but rather represents a 
general comparison of the spray cooling of surfaces using water containing dissolved 
gases against degassed water. In the next chapter a model for sparse spray cooling is 
presented and is based on degassed water. 
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Figure 17 - Droplet Dispenser 
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Figure 18 - Comparison of Water Droplets 
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Figure 19 - Transient Temperature Plot: T,, = 111 °C 
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Figure 20 - Transient Temperature Plot: T,, = 131 °C 
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Figure 21 - Transient Temperature Plot: T., = 151 °C 
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Figure 23 - Transient Temperature Plot: T,, = 182 °C 
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Figure 24 - Contour Plot: T.,= 111 °C, Time = 30 seconds 
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Figure 25 - Contour Plot: T,,= 111 °C, Time = 300 seconds 
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Figure 26 - Contour Plot: T,,= 111 °C, Time = 600 seconds 
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Figure 27 - Contour Plot: T,,= 151 °C, Time = 30 seconds 
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Figure 28 - Contour Plot: T,,= 151 °C, Time = 300 seconds 
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Figure 29 - Contour Plot: T,,= 151 °C, Time = 600 seconds 
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Figure 30 - Contour Plot: T,,= 162 °C, Time = 30 seconds 
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Figure 31 - Contour Plot: T,,= 162 °C, Time = 300 seconds 


74 


meters 


Dissolved Gases: G = 0.97 gms, Average Temperature = 126.74 °c 


0.04 


0.03 


0.02 


0.015. 0.02. 0.025 0.03 0.035 0.04 0.045 0.05 
meters 


Degassed: G = 0.97 gim2s, Average Temperature = 122.56 °c 
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Figure 32 - Contour Plot: T,,= 162 °C, Time = 600 seconds 
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Table 7 - Raw Temperature Data: G = 0.5 g/m’s 
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Table 8 - Raw Temperature Data: G = 0.96 g/m’s 
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Table 10 - Raw Temperature Data: G = 0.97 g/m’s 
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Table 11 - Raw Temperature Data: G= 1.5 g/m’s 
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Table 12 - Raw Temperature Data: G = 1.6 g/m’s 
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6. SPARSE SPRAY COOLING MODEL 

The spray cooling of a hot surface is used in many engineering applications. 
Spray cooling involves the evaporation of water droplets and, due to the high latent 
heat, represents a way to remove large amounts of heat from a surface. 

In the field of fire suppression, spray cooling from sprinklers is often used for 
extinguishment purposes as well as for the protection of surfaces. A model for the 
sparse spray cooling of a radiantly heated, semi-infinite, low thermal-conductivity solid 
surface is now presented. This model provides a relatively simple means of 
determining the transient thermal behavior of a surface that is being cooled by a sparse 
spray of water. Specifically, the code calculates transient average surface temperatures 
and spatial surface temperature distributions. The surfaces under consideration here 
are below the Leidenfrost temperature. 


6.1 Development of the Computer Model 

A previous computer model for the sparse spray cooling of a surface was 
developed by Liao [7]. This computer model predicted transient average surface 
temperatures and spatial temperature distributions for surfaces that are heated by 
conduction from below. Both a high thermal-conductivity surface (aluminum) and a 
low thermal-conductivity surface (macor) were modelled. 

The computer model presented here differs from the previous one in several 
ways. The main features that distinguish the current code from the previous code are 
listed below. 
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(1) model is for the evaporative cooling of a surface that is radiantly heated 
from above 


(2) includes accurate description of droplet distribution on the surface 


(3) closed-form solution replaces database to provide temperature drops for 
locations within 5 droplet radii from a droplet 


(4) validation of the model performed for the specific experimental setup of 
Dawson [9] 


A listing of the computer code, entitled PROGRAM MDEC, can be found in the 
appendix . 
6.1.1 Assumptions and Governing Equations 

This section contains the overall equations that govern the cooling of the 
surface and describe the physical problem. Several assumptions are made that lead to 


the final form of these equations. The governing equation and boundary and initial 


conditions for the semi-infinite solid are 


oye 
Bp ts T (16) 
T= Tyo" 3o t=0 (17) 
Ss 
-k OT g z=0,0<r<R (18) 
Sdz Cc J 
tO Th ae Z 
k= Ge aN) BY a I>R, 2=0 (19) 
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-k oT _g Z7-o (20) 


== rr, Z<0 (21) 


The initial condition of each droplet is the surface temperature corresponding 


to the location and time at which a droplet lands. This condition can be written as 


tae (22) 


jso 


@ location j,t=t j=1,ND 

With T;,, known, the temperature distribution associated with droplet j can be 
determined throughout its evaporation as was shown in Chapter 4. Therefore, for each 
droplet that impinges the surface the corresponding surface temperature prior to 
impingement must be calculated. The assumption required by this condition is that the 
temperature is uniform in a small region where the droplet is deposited. For a small 
droplet size and small temperature gradients, the assumption is valid. However, near 
the edge of a droplet the deMipeeiiis profiles shown in Chapter 4 depict large 
temperature gradients. For droplets landing near other droplets, the overestimation and 
underestimation of the droplet's initial surface temperature is expected to average out 
and cause no significant errors in predicting the overall surface cooling behavior. 

In the computer model, the quantity of interest is a temperature difference on 
the solid surface. In general, equations (16) through (21) are written in terms of the 


following quantity 
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Zz 
U(r, 2, t) =T(r, 2, t) ~Tyot 22 


s 


(23) 


where on the surface (z = 0), U is simply the difference between the temperature at a 


location and time and the initial surface temperature. 


Equations (16) through (21) take the following form when written in terms of 


6.1.2 Superposition 
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(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


The governing equations and boundary and initial conditions developed in the 


previous section describe the evaporation of a single droplet on a surface. In the 
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sparse spray cooling of a surface, several droplets will be evaporating, and thus 
removing heat, on the surface. To determine the temperature drop at a location on the 
surface, the effects of all the droplets on the surface must be superimposed. The 
governing equations and boundary conditions can be written in terms of the effects of 


| each droplet. First, the temperature difference, U;, is defined as 


ere 
Oe a ae (30) 


where T,,,, is the initial, uniform surface temperature at the deposition location of 
droplet j and T; is the surface temperature at a given location due to the effect of 
droplet j. The governing equation for the solid and the initial and boundary conditions 
given in Equations (24) through (29) can then be written in terms of U;. 
Superposition is used to obtain the total temperature difference, U, due to the 

contributions from all evaporating drops. In this way, 

oe u, (31) 
Substituting equation (31) into equations (24) through (39), the surface temperature at 


any point on the surface and at any time will have the following form 


ND 
Py Key aT eee Oy) ae eT) (32) 
j=l 


Two core assumptions are made when employing superposition. The first, 
which has already been stated, is the assumption of a uniform temperature underneath 


the droplet. The second assumption is that the evaporation of a droplet depends only 
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on the initial temperature which is associated with the droplet. That is, droplets 
landing on the surface after a droplet, j, has landed have no effect on droplet j, while 
droplets landing before only affect the initial surface temperature, T;,,, which governs 


the droplet's evaporation. 


6.1.3 Calculation Procedure 

A simplified flow chart of the computer code is displayed in Figure 33. The 
code calculates the temperature distribution and the average temperature for a portion | 
of the entire region on which droplets are impinging. The calculation domain is 
divided into a grid and a temperature is determined at each grid point. This 
temperature represents the effects due to all droplets which have landed on the surface 
up until a particular time. It is important to realize that the effects of droplets that 
have fallen outside of the calculation domain are also accounted for in the temperature 
calculation. The dimensions of the ma are such that the spatial resolution matches 
that used to obtain the experimental data, that is 0.00013 meters/pixel. 

At any time during the cooling of the surface, the temperature distributions 
and average temperature can be determined. At a grid point, the distance to each 
droplet that has impinged on the surface is calculated. Droplets within a distance of 5 
droplet radii from the grid point are considered to be in the near-field. The 
temperature drops corresponding to droplets in the near-field are calculated using the 
closed-form solution presented in Chapter 4. Thus, the temperature drop has the form 
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with the term (t - T) is set equal to 0 for t less than or equal to T. 

Droplets farther than a distance of 5 radii from the grid point are considered to 
be in the far field. The temperature drop associated with droplets in the far field is 
obtained using an instantaneous point sink solution developed by Carslaw and Jaeger 


[14]. The solution has the following form 


* 


foe meaed sitar wel ai (34) 
8 (mat) 3/2 7 
The term Q* is the strength of the instantaneous point sink. For the case of a radiantly 


heated surface, the strength is given as 


or--28e (35) 
P>Cs 


The factor of 2 in Equation (35) accounts for the fact that the droplet acts in a semi- 
infinite space rather than an infinite space. Q. is the total conductive heat that leaves 
the surface and acts to vaporize the droplet. 

The solution given by (34) assumes an instantaneous point sink ,the effects of 
which are generated at a time, t. However, an evaporating droplet removes heat over 
its entire evaporation time. For this reason, the time, t, is modified by determining 
when the average heat removal occurs during the droplet's evaporation. Figure 34 


shows for various initial surface temperatures, the average conductive heat removal as 
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a function of normalized time. The curves were asa using PROGRAM EVAP 
[10]. Referring to Figure 34, the average conductive heat removal over the 
evaporation time occurs at a time equal to 60% of the total evaporation time. 
Therefore, a delay in the deposition time of 0.6T occurs for the instantaneous point 


sink. Equation (34) is rewritten as 


Fponolavab aQtlog atte wi There (36) 
8(nmae('t-026%)) 3/7 


When using the far-field solution, any droplet that has not been on the surface 
for longer that 60% of its evaporation time has a Zero temperature drop associated 
with it. 

The choice to use the far-field solution at distances greater than 5 radii is made 
based on two factors. At approximately 5 droplet radii, the near-field solution exhibits 
some computational oscillations for times prior to the droplet's evaporation. The 
oscillations result from the oscillatory nature of the integrand and the numerical 
solution to the integral. Though the resulting oscillations are found to be less than 
0.15 ° C, they are nonetheless undesirable. Additionally, a radius of influence of 
approximately 3 to 4 droplet radii has been determined form previous studies [16] and 
was discussed in Chapter 5. This suggests that at a distance of 5 radii from a droplet, 
the cooling effects are felt only in a remote sense and the far-field solution provides 
this type of result. Figure 35 shows the difference between the results obtained using 
the far-field solution and the near-field solution at a distance of 5 radii for various 


initial surface temperatures. The figure shows that the difference is always contained 
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in a narrow band of 0.15 °C. 


6.1.4 Droplet Distribution 

Prior to calculating the temperature at the grid points, information on the 
droplets that have impinged the surface over a specified time period is obtained using 
a droplet distribution subroutine (SUB DDIST). This subroutine provides the location 
of the droplet on the surface, the deposition time, the temperature at the point where 
the droplet landed, and the evaporation time of the droplet. Of these pieces of 
information, the positions of the droplet on the surface needed clarification with 
respect to the initially proposed random distribution. 

In the sparse spray cooling experiments [9], attempts were made to distribute 
the droplets in a random fashion over a circular area of the solid surface. Such a 
random distribution was thought to simulate a spray cooling process more accurately. 
To this end, a positioning mechanism that consists of three solenoid-controlled 
bumpers was built to collide with and impart motion to the vertically-suspended 
droplet generator. The droplet generator, which hung form wires, would swing in a 
predetermined area as it impacted with the various bumpers. To keep the motion from 
decaying or falling into a particular pattern, a motorized cam was used to periodically 
pluck one of the suspension et 

In the previous sparse spray cooling model, Liao 7] attempted to simulate the 
droplet distribution by assuming a constant velocity distributor moving within a 


circular area. Upon hitting the perimeter, the distributor bounces back on a new 
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trajectory according to a perfectly elastic collision. 

Keeping in mind the ideal random droplet distribution, initial versions of the 
current model used a random number generator to assign the droplet positions. 
Positions were assigned in polar coordinates with a random angular and radial 
position. To ensure an evenly-distributed droplet density over the entire circular area, 
the square root of the random number was taken and assigned as the radial location. | 

After preliminary testing of the model, an attempt was made to more accurately 
describe the droplet distribution that occurred in the laboratory. Figure 36 shows the 
droplet distribution recorded during a typical experiment. The distribution 
corresponds to the area that is viewed by the infrared camera, with the outside radius 
equal to 3 cm. The total area that droplets impinge upon is approximately 3.25 cm in 
radius. Considering a circular area with a radius of 2 cm within the distribution 
shown in Figure 36, the ratio of the number of droplets that fall within this area to 
the total number of droplets shown is 74 to 1. Thus, the assumption of a uniform 
droplet density over the wetted area is incorrect. Clearly, the density is higher for the 
inner circular area, which represents about half of the total area. 

In order to determine the experimental distribution, the first step was to 
characterize the motion of the droplet generator. As mentioned previously, the 
generator swings back and forth in a region that is confined by a set of bumpers. The 
bumpers move in and out in the horizontal plane in a synchronous fashion. During 
their motion, the bumpers have an outermost position that corresponds to a radius of 


3.25 cm and an innermost position that corresponds to a radius of roughly 1.8 cm. 
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These positions suggest that the droplet generator can move freely only in a 1.8 cm 
radius. The generator is not expected to have a radial position of the maximum 3.25 
cm because the bumpers will impede the generator as they move back toward the 
center after reaching this radius. Studying the distribution shown in Figure 36 
confirms the generator motion expected from the experimental setup. 

To model the droplet distribution, a relationship between the random number 
produced by the random number generator and the radial location needs to be 
determined. In this way, the randomness of the droplet position is maintained while 
distribution characteristics are accounted for. diMarzo [17] realized that the integral of 
the function describing the droplet distribution constituted the aforementioned 
relationship. The function describing the droplet distribution has the following 
boundary conditions based on a normalized radius. 
d=0 @r=0: 
ae a iors the correct random distribution for the free area of motion 


d=0 @r=€ where € is the radius of free motion, requires a continuous function 
with a maximum at this radius 


d=0 @r=1: requires generator to never reach outermost position 
J,! d dr = 1: requires all droplets to fall in the given area (normalized function) 
The form of the function d, then, is 


Ce le Tie iO T, seat ACT cent any thie (37) 


With € = 0.56 (~ 1.8/3.25), the coefficients in (37) can be solved for to obtain 


Integrating (38) yields the desired relationship between the random number, D, and the 


a3 


dqvavo4iort "a" O7VGArey + 11.497* ee or (38) 
radial position, r, of a droplet. Therefore, 
Dre= le 8327 Faso teSu.6 64, tee 8 8S Lot 22 (39) 


Note that D(0.707) = 0.76. This quantity represents the ratio of the number | 
of droplets falling in an area equal to half that of the entire impingement area to the 
total number of drops. The ratio obtained using Figure 36 is 0.74, indicating that the | 
function d given by (38) provides a reasonable representation of the experimental 
droplet distribution. To further check the validity of the function given in Equation 
(39), the percentage of droplets falling in any particular radius in Figure 36 can be 
calculated. By dividing the area of droplet impingement into evenly-spaced concentric 
circles, a range of percentages corresponding to each radial location can be calculated. 
The ranges are determined based on two positions of the outermost concentric ring 
with respect to the outermost droplet. These two positions are determined such that 
the outermost droplet lies along either the inner radius or the outer radius of the last 
concentric ring. Figure 37 shows the relationship between the random number, D, and | 
the normalized radial position for the distributions given a linear relationship, a spatial 
relationship, and the relationship given in (39). The droplet percentage range at 
incremental radial locations is also plotted. The fact that the function given in 
Equation (39) falls within these ranges further supports that the droplet distribution has 
been characterized correctly. A typical droplet distribution computed by the 


subroutine SUB DDIST is shown in Figure 38. The distribution is for a mass flux of 
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0.5 g/m’s. 


6.2 Results and Analysis — 

The results of the sparse spray cooling model and sparse spray cooling 
experiments [9] are found in Figures 39 through 55 , with raw data contained in 
Tables 13 through 16. Dawson's experimental data includes initial surface 
temperatures ranging from 110 °C to 182 °C and water mass fluxes ranging from 0.2 
g/m’s to 1.7 g/m’s. Calculated results are compared against four data sets 
corresponding to the degassed water experiments. The cases selected for comparison 
represent a relatively wide span of both initial surface temperature and mass flux that 
are found in evaporative cooling. Large experimental mass fluxes that produced 
steady state temperatures below 100 °C were not used for comparison because these 
temperatures are not relevant to evaporative cooling applications. Additionally, 
experiments with initial surface temperatures of 182 °C could not be used because 
nucleate boiling occurs at this temperature. The four experiments selected allow for a 
comparison of both mass flux and initial surface temperature effects. Thus, the results 
are for initial surface temperatures of 131 °C and 162 °C, each having a mass flux of 
0.5 g/m’s, and for initial surface temperatures of 151 °C and 162 °C, having mass 
fluxes of 0.96 g/m’s and 0.97 g/m’s, respectively. 

Although only four data sets are selected to validate the code, it is important to 
realize how a data point is obtained. Each infrared snapshot of the surface represents 


a grid of approximately 354 pixels x 262 pixels. A data point is calculated by 
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averaging every fifth pixel intensity values, or approximately 3710 pixel values. 
Transient average surface temperatures are comprised of data points at every 50 
seconds for 24 minutes. Therefore, each data set contains averages of approximately 


107,000 pixel values. 


6.2.1 Transient Average Surface Temperature 

In Figure 39 through Figure 42 the average surface temperature is plotted 
versus time for the four cases listed in Section 6.2. Experimental data is plotted using 
open symbols while the calculated results are plotted using closed symbols. With the 
exception of the results corresponding to an initial surface temperature of 162 °C anda 
mass flux of 0.5 g/m’s, the calculated and experimental results agree well. The code 
is able to predict steady state temperatures that agree with those found in the 
experiment. Additionally, the code is able to capture the transient behavior of the 
surface. 

Experimental data for an initial surface temperature of 151 °C and a mass flux 
of 0.96 g/m’s does not exist past a time of approximately 15 minutes. The code was 
run for a longer time, however, in order to observe the predicted steady state behavior. 
The initial portion of the calculated results agrees with the experimental data and it 
appears that if experimental data existed for longer times it would follow the trend 
observed in the calculated results. The largest disagreement between the calculated 
results and the experimental data is found for the case having an initial surface 
temperature of 162 °C and a mass flux of 0.5 g/m’s. In the first five minutes of 
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cooling, the experimental temperatures are lower than the predicted temperatures. At a 
time of five minutes, the experimental data exhibits a discontinuity where the surface 
seri have overcooled, recovered slightly, and cooled again. At ten minutes into 
the cooling, the surface appears to begin to approach the steady state temperature 
predicted by the code. A simple linear extrapolation of the experimental data points 
indicates that the surface would reach its steady state temperature at toni of 
approximately 58 minutes. The most likely reason for the deviant behavior of this 
data set is that an initial surface temperature of 162 °C is near the range of nucleate 
boiling. The data suggests that nucleate boiling may have occurred during the initial 
cooling period causing the surface to cool to lower temperatures than expected. After 
the initial sharp cooling of the surface, it tm that the continual impingement of 
water on the surface causes a suppression of the nucleate boiling. The suppression 
allows for a slight recovery of the surface and for the remaining of the surface cooling 
to be governed by evaporation. The deviant behavior is not seen in the data | 
corresponding to an initial surface temperature of 162 °C and mass flux of 0.97 g/m’s 
because the relatively high mass flux is able to suppress the nucleate boiling. 

From the experimental results shown in Chapter 5 and contained in reference 
[9], two important parameters that govern the cooling of the surface are evident. 
These parameters are the steady state surface temperature and the behavior of the 
surface as it wre to the steady state temperature. The latter parameter is described 
by a time constant. Various quantities, such as the mass flux, the initial surface 


temperature, or the solid material properties, may affect these two parameters. It is 
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important to determine any relationships between these parameters and the quantities 
that may affect them. 
A maximum temperature difference is defined in the following way 


AT. =T..-Tss (40) 


Dawson [9] found a relationship between the surface temperature and time of the 


following form 


_£ 
T= AT,€ *+T,, (ag) 


A normalized temperature can now be defined as 


_ T-Tss 
Ss AT, . (42) 


In terms of ©, equation (41) becomes 
@ = e7t/t (43) 


In this equation, T, is the time constant mentioned previously. 

The experimental and calculated data are plotted in the form of © versus t in 
Figure 43. The steady state temperature used to plot the data is obtained using the 
calculated results. With the exception of the experiment consisting of an initial 
surface temperature of 162 °C and mass flux of 0.5 g/m’s, the experimental and 
calculated data collapses. The transient results plotted in this form indicate that the 
transient thermal behavior of the surface is independent of mass flux and initial 
surface temperature. The calculated results corresponding to an initial surface 
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temperature of 162 °C and mass flux of 0.5 g/m’s are consistent with the other sets of 
data, indicating that the experimental data for that case represents an anomaly and that 
the experimental conditions may not have been conducive to evaporative cooling. 
Figure 43 gives both theoretical and experimental support to the conclusion that the 

| time constant, T, is independent of mass flux and initial surface temperature. 

To determine the time constant, a curve is fit to the calculated results in the 
form of e*. To do this, In® is plotted against t and a linear regression is performed. 
A regression is performed only for the initial transient (t < 3.5 min) because this is the 
region where calculated results have the least scatter. A value of T equal to 230 
seconds was calculated in this way. 

Since T was found to be independent of the mass flux and initial surface 
temperature, one might expect that it instead depends on the properties of the solid. 
At steady state, cooling effects are expected to be felt within a certain depth of the 


solid. This depth is called the penetration depth, 1, and is given as 
T (44) 


Using o, = 5.79*10”’ m/s for macor and tT = 230 seconds as calculated, | is found to 
be 1.15 cm. A radius of influence was discussed in Section 5.2.1 and determined to 
be between 1.05 cm and 1.5 cm for the experimental conditions under consideration. 
Therefore, the model supports the experimental data in suggesting that the time 
constant is a function of the solid properties, and that the penetration depth defined by 
(44) is of the same order of magnitude as radius of influence of a single droplet on the 
surface. 
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6.2.2 Surface Temperature Contours 

A more qualitative validation of the sparse spray cooling code is found in the 
temperature contour plots contained in Figures 44 through 55. The contour plots are 
simply snapshots of the surface taken at three different times during each experimental 
run. The times are chosen to represent the early portion of the surface cooling (t = 50 
seconds), an intermediate portion of the surface cooling (t = 300 seconds), and a 
steady state portion of the surface cooling (t = 600 seconds). Again, the cases used 
for the qualitative validation are those described in the previous section. Two contour 
plots, corresponding to experimental results and calculated results, respectively, are 
shown in each figure. 

Several features are similar in both the experimental and calculated contour 
plots. The first feature is that both plots contain local concentric rings of contour lines 
indicating locations of evaporating droplets. At a given time, the number of droplets 
indicated by the contours are approximately the same for both experimental and 
calculated results. There may be slight differences in the number of droplets due to 
the fact that the snapshot represents only a portion of the surface and droplets may fall ; 
outside of the displayed region. Another observable trend seen in both sets of results 
is that the number of droplets increases with time. This is expected because the 
evaporation times are longer as the surface cools. 

The temperatures indicated on the contours also agree between the experiments 
and the code. The innermost contours have the lowest temperatures while the 
outermost have the highest. The innermost contour temperatures for the experiments 
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are generally lower than those for the code. However, these contours represent 
regions underneath the droplets where the infrared camera yields erroneous data. The 
more important temperatures are those associated with the contours that run through 
the main portion of the plots. These contours indicate how the surface behaves as a 
whole, not just in the local region where droplets exist. These temperatures are 
consistent between the experimental and calculated results. A decrease in these 
temperatures occurs as time increases and the surface cools. 

The main difference between the experimental and calculated contour plots is 
that the calculated contours appear smoother than the experimental contours. The 
jagged contour lines in the experimental results occur because of the fluctuation in the 


intensities at the various surface points. 


6.3 Summary 

A model for the prediction of the transient thermal behavior of a low thermal 
conductivity, semi-infinite solid subjected to a random array of water droplets has 
been presented and tested against experiments. The model uses superposition and 
closed-form solutions to calculate surface temperature drops caused by all droplets that 
have impinged the surface. A detailed description of the droplet distribution that 
occurred during experiments is incorporated in the model. | 

Validation of the model against experimental data shows that the model is able 
to predict the transient thermal behavior of the solid both quantitatively and 
qualitatively. The code is able to calculate the steady state temperature of the surface 
given an initial surface temperature and a mass flux and is also able to predict how the 
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surface reaches its steady state temperature. Furthermore, it was found that a 
dimensionless temperature could be defined such that time plots of both the calculated 
and experimental data in this form results in a collapsing of all the data. From this 
fact, the time constant, which describes the transient cooling of the surface, is 
determined to be of the mass flux and is found instead to depend on the properties of 
the solid. 

Perhaps the greatest source of error expected in the code is the use of the 
closed-form solution in the near-field. Figures 15 and 16 in Chapter 4 show that 
relatively large discrepancies in temperature occur between the closed-form solution 
and the single droplet model. To evaluate the error caused by these discrepancies, an 
average difference in temperature of 10 °C is assumed. Underneath the droplet 
differences are usually from 5 °C to 15 °C. The contour plots in show that at steady 
state there are generally between 5 and 8 droplets existing on the surface. The 


difference expected in average temperatures is given by 


| (45) 
peg Bezida) 


ave A, 
where A, is the area underneath the a droplet (2.77*10° m”), A, is the total area of 
calculation (1.5*10° my), and AT is approximately 10 °C. Using (45) ATuye is 
calculated to be 0.185 °C. Since the experimental infrared thermography itself 
produces errors that are +2 gray shade and thus +2 °C, the error in using the closed- 
form solution is negligible in the calculation of the average surface temperature. 
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Figure 33 - Simplified Program Flow Chart 
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Figure 34 - Normalized Average Conductive Heat Removal 
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Figure 35 - Near and Far Field Solution Difference 
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Figure 36 - Typical Experimental Droplet Distribution 
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Figure 37 - Relationships Between D and Radial Position 
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Figure 38 - Typical Calculated Droplet Distribution 
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Figure 39 - Transient Temperature Plot: T,, = 131 °C 
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Figure 40 - Transient Temperature Plot: T., = 162 °C 
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Figure 41 - Transient Temperature Plot: T,, = 151 °C 
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Figure 42 - Transient Temperature Plot: = 162 °C 
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Figure 43 - Experimental and Calculated Transient Dimensionless 
Temperature Plot 
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Experiment: G =0.5 g/m or Average Temperature = 125.12°C 
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Calculated: G=0.5 g/m?s, Average Temperature = 125.26°C 
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Figure 44 - Contour Plot: T,, = 131 °C, Time = 50 seconds 
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Experiment: G=0.5 g/m s. Average Temperature = 114.25 °C 
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Calculated: G = 0.5 g/m2s, Average Temperature = 114.99 °C 
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Figure 45 - Contour Plot: T,, = 131 °C, Time = 300 seconds 
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Figure 46 - Contour Plot: T., = 131 °C, Time = 600 seconds 
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Experiment: G=0.5 g/m’ s, Average Temperature = 154.61 °C 
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Figure 47 - Contour Plot: T,, = 162 °C, Time = 50 seconds 
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Experiment: G=0.5 g/m 2, Average Temperature = 148.10 G 
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Calculated: G=0.5 g/m 2. Average Temperature = 151.45 °C 
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Figure 48 - Contour Plot: T,, = 162 °C, Time = 300 seconds 
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Experiment: G = 0.5 g/m?s, Average Temperature = 138.93 °C 
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Calculated: G=0.5 g/m 2, Average Temperature = 147.06 ae 
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Figure 49 - Contour Plot: T,, = 162 °C, Time = 600 seconds 
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Figure 50 - Contour Plot: T,, = 151 °C, Time = 50 seconds 
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Experiment: G = 0.96 g/m 2, Average Temperature = 117.36 C 
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Calculated: G= 0.96 gim2s, Average Temperature = 121.2 °c 
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Figure 51 - Contour Plot: T,, = 151 °C, Time = 300 seconds 
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Experiment: G = 0.96 gim’s, Average Temperature = 119.49 °C 
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Calculated: G= 0.96 g/m2s, Average Temperature = 111.55 °.C 
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Figure 52 - Contour Plot: T,, = 151 °C, Time = 600 seconds 
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Experiment: G = 0.97 g/m s, Average Temperature = 151.38 °C 
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Calculated: G = 0.97 g/m2s, Average Temperature = 153.56 °C 
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Figure 53 - Contour Plot: T,, = 162 °C, Time = 50 seconds 
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Calculated: G = 0.97 gms, Average Temperature = 135.96 °C 
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Figure 54 - Contour Plot: T,, = 162 °C, Time = 300 seconds 
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Experiment: G = 0.97 oNaes. Average Temperature = 122.56 7C 
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Calculated: G = 0.97 g/m2s, Average Temperature = 125.90 °C 
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Figure 55 - Contour Plot: T,, = 162 °C, Time = 600 seconds 
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Table 13 - Raw Temperature Data: G = 0.5 g/m’s 
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Table 14 - Raw Temperature Data: G = 0.5 g/m’s 
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Table 15 - Raw Temperature Data: G = 0.96 o/m’s 


128 


| aed Time | [ Experimen ime alculated 
nit) | © | mm | | 
[0.00] 16200] 0.00 | UU 
e+e et ee 
eee fe + ee 
p 150 | 143.440 | 250 [14453 | 
ageib0 del topmnls9.0 few waht 4 aphsiees 190 Orie] 
p3.00) | 133.52) | 5.00 | 135.98 
P 3.50 | 137.56 5.83 135.69 
[es0UT te wlooloes | O.Ole |e ot oOn 
is: Supe ain 34 22m] n7-50 ims |e SOT | 
| S00 a | PZ 7.32 B.33-" [PI 28.96. 
[ann 0s [tune 32.09. 0esn an 0.A7 wan [aan 109.07 Go| 
ee 
ig ges |e od Aanoe dT | ean 2649 cal 
Ss ol a eid iam A ees Dial 2 Biel 
(pee 00 ame [pee 24 SON mS. GS] ag 125,.S80 en | 
al el SN le 
, 9.00 | 12065 | 15.00 { 12335 | 
STO ase | eT tiase | 
Se ee 
moon paniz0eses [ess Tigtt = | 
jimi, Opa | peal 22 Coupes | and 01 7am | wand 1 7.86 | 
[12.00 | 120.84 | 2000 | 11627 | 
Pr 2 a Be 
| 13.00 [| 11860 | 21.67 | 11452 | 
| 13.50 | 12272) | 2250 | 11655 
| 14.00 [| 123.84 | 23.33 [ 11609 | 
TIE Ee aed ed ee 
| 15.50 | 12036) | 25.83 [| 11266 
PEE Re De 
| 1650 | 12546 | 27.50 | 11480 | 
| 17.00 | 12243 | 28.33) | 110.02 
BA as Ee ee ee eee eee 
BAB OC at [120.55 jis | eee ee eae 
SOU Me 121.83. iis fet Si tape Ds, 
ees SOD RE eet cO. 440 eer] = men near wereree eee ee 
2710.50 ee WRI20.83e5s | ole ea 
20-00) often 24-0tw biwengob 20 [ae Ones & 
PE DEP CE en ea 
621.00 Be ead ISS firiy nhttew wey lainey alonio 
Perce el ae 
22.00) See ipa 20.42 Say ey ee | ea | 
eeu ee eee iee VON ee S| ee 
PEEL Coy A ee 
Table 16 - Raw Temperature Data: G = 0.97 g/m’s 
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7, CONCLUSIONS 

The development of the transient droplet geometry that is used in a model for 
the prediction of a single droplet evaporating on a low thermal-conductivity, semi- 
infinite solid subjected to radiant heat input was presented. Using the droplet 
configuration model, it was found that while the shape factor has a strong effect on 4 
evaporation time, the initial solid-liquid-vapor contact angle has virtually none. 

With the model for the droplet geometry incorporated and a modification to the. 
radiation terms to include the geometric effects of the droplet, an in-depth validation 
of the single droplet model against experiment was performed. - Validation of the 
model showed that it is able to accurately predict transient surface temperature profiles 
caused by the evaporation of a droplet. 

Experimental work using water containing dissolved gases was performed. The 
results were compared against previous experiments using degassed water. In general, 
there appeared to be no major differences in the transient cooling of the surface or the 
steady state temperature reached by the surface. Qualitative results showed 
differences, specifically with regards to the size of the wetted region undemeath the 
droplet and temperature gradients near the droplet's edge. 

A model for the prediction of surface temperatures of a hot, semi-infinite solid 
being cooled by a sparse spray of degassed water has also been presented. The 
validation of the single droplet evaporation model that includes the detailed 
descriptions of the transient droplet geometry and of the radiation effects led to the use 


of a simplified, closed-form solution in the sparse spray cooling model. This solution 
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decreases run times and memory usage by eliminating the need for a database. 

The model is able to predict the transient thermal behavior of the solid surface 
over a range of initial surface temperatures. The surface temperatures of interest for 
fire protection applications are generally between 100 °C and 160 °C. 

Using results obtained both experimentally and through the code, it was found 
that when plotting a dimensionless eee versus time, all of the data collapses 
onto a single curve. In this way, the independence of the time constant that governs 
the transient cooling of the surface and the mass flux is established. Results support 
Dawson's hypothesis [9] that the properties of the solid are the factors which 
determine the time constant. Using the time constant to calculate the penetration 
depth yields a value equal to the radius of influence of a single droplet evaporating on 
the surface. | 


Aleit 


APPENDIX - PROGRAM MDEC 
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... SUSANNE C. TINKER 

.. MASTER'S THESIS - 1992 - 94 

...ADVISOR: DIMARZO 

.. UNIVERSITY OF MARYLAND 

..PROGRAM TO CALCULATE THE TRANSIENT SURFACE TEMPERATURE DISTRIBUTION 
«CAUSED BY SPRAY COOLING OF A SURFACE 


COMMON TDROP, RADO, TSO 


REAL LTAU 


REAL SAREA, V0, BETAO, RHOW, RHOS, CS, DELTA, TCP, ALPHAS, LAMDA, KS, PI 
REAL LDX, LDY, RDX, RDY 

REAL XWIDTH, YWIDTH 

REAL U, UTOT, DIST, TIME, TBEGIN, STEP, UCRIT 

INTEGER NX, NY, NDATIM, L J, K, M, KLAST, LNDROP 

REAL FREQ, DELTIME, G, TEND; TTOT, TAVE 

INTEGER TDROP 

REAL RADO, TSO 


DIMENSION ARRAYS 


REAL XCALC(100), YCALC(100), T(100,100), DATTIM(50) 


.... DIMENSION ARRAYS CONTAINING DROPLET INFORMATION 


INTEGER NDROP(1000) 
REAL XD(1000), YD(1000), DEPTIME(1000), T0(1000), TAU(1000) 
CHARACTER*1 SCFLAG(1000) 


... THE AREA OF THE WETTED SURFACE AREA (TESTED AGAINST EXPERIMENT .0033 M‘2) 


PARAMETER (SAREA = 0.0033) 


.... THE INITIAL SHAPE PARAMETER BETAO (TESTED AGAINST EXPERIMENT FOR 2.3) 


PARAMETER (BETAO = 2.3) 


... DEFINE THE CALCULATION DOMAIN (ONLY A PORTION OF THE WETTED AREA) 


LDX = 0.0105 
LDY = 0.0155 
RDX = 0.0545 
RDY = 0.0495 
XWIDTH = RDX - LDX 
YWIDTH = RDY - LDY 


..DETERMINE THE MESH SIZE TO BE USED FOR THE CALCULATION DOMAIN 
PRINT *, 'ENTER THE NUMBER OF GRID POINTS IN THE X-DIRECTION' 
..READ (6,*) NX 

...PRINT *, ENTER THE NUMBER OF GRID POINTS IN THE Y-DIRECTION’ 


C..... MORE VARIABLES TO BE INPUT BY USER THROUGH KEYBOARD 
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C.....PRINT *, ENTER THE MASS FLUX (IN KG/M’2-S) ' 
C.....READ (6,*) G 
G = 0.00097 
C.....PRINT *, ‘ENTER THE INITIAL SURFACE TEMPERATURE (IN CELSIUS) ' 


C.... PRINT *, ENTER THE NUMBER OF TIMES TO OUTPUT DATA (INTEGER 50 max) 
C.....READ (6,*) NDATTM 
NDATIM = 13 
C.....PRINT *, "ENTER THE TIMES TO OUTPUT THE SPATIAL DATA (50 MAX) 
C.....DO 51=1, NDATIM 
C.....READ (6,*) DATTIM() 
C..... DATTIM(1) = 100. 
C.....PRINT *, ‘ENTER TOTAL CALCULATION TIME" 


RADO = BETAO * (3.0 * VO¥(4*3.141592654))**(1,3.) 


C..... CALCULATE THE TIME STEP AND THE FREQUENCY OF THE DROPLETS 
FREQ = (G * SAREA) / (VO * RHOW) 
DELTIME = 1/FREQ 
TDROP = INT(FREQ*TEND) 


CALL DDIST(TEND, DELTIME) 


e 

C 

C.....OPEN FILES 
OPEN(11, FILE = 'dropdat.dat', STATUS = 'OLD’) 
OPEN(12, FILE = 'tempdat.dat’, STATUS = NEW) 
OPEN(13, FILE = 'y.dat’, STATUS = 'NEW) 
OPEN(14, FILE = 'x.dat’, STATUS = ‘NEW) 
OPEN(15, FILE = ‘tempave.dat’, STATUS = "NEW) 


C.....READ IN DROPLET DATA FROM DROPDAT.DAT AND INITIALIZE FLAG TO 
C.....STOP CALCULATION IF TEMPERATURE OF CENTER OF DROPLET IS 


DO 40 K=1, TDROP 
READY RE NDROP(K), DEPTIME(K), TOK), TAU(K), XD(K), YD(K) 


0 
40 CONTINUE 
C.....CALCULATION OF TEMPERATURES AND AVERAGE SURFACE TEMPERATURES AT VARIOUS TIMES 


DO 80 TIME = 50., TEND, 50. 
WRITE(12,*) TIME 
IF(TIME .LT. TEND) THEN 
LNDROP = INT(FREQ*TIME) 
ELSE 
LNDROP = INT(FREQ*TEND) 
END IF 


C..... TEST TO SEE IF DROPLET IS STILL AFFECTING SURFACE 
DO 60 M = 1, LNDROP 
DIST = 0.0 
CALL NEAR(TIME, M, DIST, DEPTIME, U, T0) 
IF (U .GT. 0.0) SCFLAG(M) ='1' 
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60 CONTINUE 


C..... INITIALIZE ARRAYS AT EACH TIME 
DO 75 J=1, NY 
YCALC(J) = 0.0 
DO 76 I= 1, NX 
XCALC() = 0.0 
TL) = 0.0 
76 CONTINUE 
75 CONTINUE 


DO 50J=1, NY 
YCALC(J) = 0.0155 + (YWIDTH/NY)*(J) 
DO 55I=1, NX 
XCALC() = 0.0105 + (XWIDTH/NX)*() 
UTOT = 0.0 
IF (LNDROP .EQ. 0) THEN 
TJ) = TSO 
ELSE 
DO 70 K = 1, LNDROP 
IF (SCFLAG(K) .EQ. 'I') GOTO 70 
IF (DEPTIME(K) .GT. TIME) GOTO 100 


DIST = SQRT((XCALC(1) - XD(K))**2 + (YCALC() - YD(K))**2) 
IF (DIST .GT. 5.0 * RADO) THEN 
CALL FAR(TIME, K, DIST, DEPTIME, U, TO) 


Cc PRINT *, 'F ', DIST, K, U 
ELSE 
CALL NEAR(TIME, K, DIST, DEPTIME, U, T0) 
C PRINT *, 'N ', DIST, K, U 
END IF 


IF (U .GT. 0.0) U=0.0 
UTOT = UTOT + U 
70 CONTINUE 
END IF 
T(LJ) = TSO + UTOT 
TTOT = T(LJ) + TTOT 
55 CONTINUE 


C..... WRITE TO CONTOUR FILES ONLY FOR SPECIFIED TIMES 


IF (TIME .EQ. 50. .OR. TIME .EQ. 300. OR. TIME .EQ. 600.) THEN 
Rated (TLD, I= 1, NX) 


WRITE(13,*) YCALC(J) 
CONTINUE 


WRITE(14,*) (XCALC(), I = 1, NX) 
100 TAVE=TTOT((NX)*(NY)) 


50 


C..... WRITE TRANSIENT AVERAGE SURFACE TEMPERATURE TO FILE 
WRITE(15,*) TIME, TAVE 
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CLOSE (UNIT = 12) 
CLOSE (UNIT = 13) 
CLOSE (UNIT = 14) 
CLOSE (UNIT = 15) 
END 


SUBROUTINE DDIST(TEND, DELTIME) 
C..... THIS SUBROUTINE CREATES A FILE CALLED DROPDAT.DAT WHICH CONTAINS ALL THE DROPLET 


C.....TO BE USED IN THE MAIN PROGRAM THIS SUBROUTINE IS ONLY CALLED ONCE DURING THE 
C.....EXECUTION OF THE MAIN PROGRAM 


COMMON TDROP, RADO, TSO 
COMMON /CONST/ ALPHAS, CS, DELTA, KS, LAMDA, PL RHOS, RHOW, TCP, VO 


REAL LTAU , RAND 


REAL X(1000), ¥(1000), DTIME(1000), TAU(1000) 

REAL TEMP(1000), TIME, DIST, THETA, G, TSO, RADO, UTOT, RADHI 
REAL TEND, DELTIME, R, U, UCRIT, RADLOC, RANDOM, RADLOW 
REAL DELTA, KS, TCP, ALPHAS, CS, V0, RHOW, RHOS, LAMDA, PI 
INTEGER NUM(1000), TDROP, K, I, J, COUNT 

CHARACTER*1 SCFLAG(1000) 


OPEN (UNIT = 11, FILE = 'dropdat.dat’, STATUS = NEW) 


C..... INITIALIZE RANDOM NUMBER GENERATOR SEED 
R=0.0 


C..... INITIALIZE COUNT 
COUNT = 0 


C..... INITIALIZE TIME 
TIME = 0.0 


C..... INITIALIZE STOP CALCULATION FLAG 
DO 151=1, TDROP 
SCFLAG() = '0 
15 CONTINUE 


DO 20 K = 1, TDROP 
TIME = TIME + DELTIME 
NUM(K) = K 


C..... DETERMINE DROPLET LOCATION ON SURFACE USING ZETA = .56 


THETA = RAND(R)*PI*359./180. 


RADHI = 1.0 
RADLOC = (RADLOW+RADHI/2. 
25 G=1.83*RADLOC**5-5.66*RADLOC**4+3.83*RADLOC**3+RADLOC**2 
IF (G - RANDOM) .GT. 0.001) THEN 
RADHI = RADLOC 
RADLOC = (RADHI + RADLOWy?2. 
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GOTO 25 
ELSEIF ((G - RANDOM) .LT. -0.001) THEN 
RADLOW = RADLOC 
RADLOC = (RADHI + RADLOW)Y/2. 
GOTO 25 
END IF 
X(K) = 0.0325*RADLOC*COS(THETA) + 0.0325 
Y(K) = 0.0325*RADLOC*SIN(THETA) + 0.0325 


C..... CALCULATE DEPOSITION TIME OF DROPLET 
DTIME(K) = TIME 


C..... CALCULATE SURFACE TEMPERATURE WHERE DROPLET HAS C.....LANDED 


IF (K EQ. 1) THEN 
TEMP(K) = TSO 
ELSEIF (K .GT. 1) THEN 


C..... CALCULATE TEMPERATURE OF SPOT AT WHICH NEW DROPLET LANDED TO DETERMINE TO FOR ALL 
C..... DROPLETS 


DO 35 J=1, K-1 
IF (SCFLAG(J) EQ. '!') THEN 
GOTO 35 
ELSE 
DIST = 0.0 
CALL NEAR(TIME, J, DIST, DTIME, U, TEMP) 
IF (U .GT. 0.0) SCFLAG(J) ='1' 
END IF 


Cc IF (SCFLAG(J) EQ. 'l’) GOTO 35 
DIST = SQRT((X(K) - X{J))**2 + (Y(K) - Y))**2) 


C.....USE FAR FIELD SOLUTION FOR DISTANCES GREATER THAN 5 RADI 
IF (DIST .GT. 5.0* RADO) THEN 


e: PRINT *, 'F ', K 
CALL FAR( TIME, J, DIST, DTIME, U, TEMP) 
C PRINT *, 'F', K, J, U 
ELSE 


C..... USE NEAR FIELD SOLUTION FOR DISTANCES LESS THAN 5 RADI 


Cc PRINT *,'"N '",K 
CALL NEAR(TIME, J, DIST, DTIME, U, TEMP) 
Cc PRINT *, ‘N', K , J, DIST/RADO, U 
END IF 


IF (U .GT. 0.0) U=0.0 
UTOT = UTOT + U 


C.....CALCULATE TIME AT WHICH DROPLET WILL EVAPORATE 
TAU(K) = LTAU(K, TEMP) + DTIME(K) 
C.....DETERMINE IF STOP CALCULATION FLAG SHOULD BE USED 


C..... WRITE TO DROPDAT FILE 


3:7 


WRITE (11,27) NUM(K), DTIME(K), TEMP(K), TAU(K), X(K), Y(K) 
27 FORMAT (1X15, 1X, F15.6, 1X, F15.6, 1X, F10.3, 1X, 2E11.3) 
20 CONTINUE 

CLOSE (UNIT = 11) 


RETURN 
END 


REAL FUNCTION LTAU(K,T) 
REAL (1000) 
INTEGER K 


C.....CURVE FIT TO EVAPORATION TIME DATA OBTAINED FROM GLENN WHITE'S CODE 
LTAU = 1300. * EXP(-0.03 * T(K)) 


RETURN 
END 


REAL FUNCTION F(X) 


COMMON /ARG/ ARG1, ARG2, ARG3 

REAL X 

REAL ARGI, ARG2, ARG3 

REAL ERF, BESIOE, BESJO, BESJ1 

F = BESJQ(ARGI * X) * BESJI(X) * (ERF(ARG2 * X}ERF(ARG3 * X)) /X 


RETURN 
END 


SUBROUTINE NEAR(TIME, K, DIST, DEPTIME, U, TEMP) 
C..... THIS SUBROUTINE CALCULATES THE TEMPERATURE DEPRESSION AT A POINT 
C.....DUE TO A DROPLET LOCATED 
C.... LESS THAN 10 DROPLET RADII AWAY FROM THE POINT 
C..... THE SOLUTION IS A CLOSED FORM SOLUTION FOR A DISK OF C.....CONSTANT STRENGTH 
C.....GIVEN BY CARSLAW AND JAEGER AND MODIFIED BY GLENN WHITE 


COMMON TDROP, RADO, TSO 
COMMON /ARG/ ARGI, ARG2, ARG3 
COMMON /CONST/ ALPHAS, CS, DELTA, KS, LAMDA, PL RHOS, RHOW, TCP, VO 


EXTERNAL F 

REAL RADO, KS, ALPHAS, QL QD, RAD, DELTA, TCP, CS, PL RHOW 
REAL RHOS, V0 

REAL TIME, U, DIST, LTIME, TEMP(1000), DEPTIME(1000), TSO, C 
REAL ARGI, ARG2, ARG3 

INTEGER K, |, J, TDROP 


C.....SPECIFY VARIABLES TO BE USED BY SUB QAGI 
REAL EPSABS, EPSREL, RESULT, ABSERR, BOUND, WORK(5000), F 


REAL LTAU 
INTEGER INF, IER, LIMIT, LENW, TWORK(1000), NEVAL, LAST 
C.....DEFINE VARIABLES NEEDED BY SUB QAGI (CMLIB LIBRARY) USED FOR CALCULATION OF SI 


INTEGRAL 
BOUND = .00001 
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C..... DEFINE LOCAL EXISTENCE TIME OF DROPLET TO BE USED IN CLOSED FORM SOLUTION 
LTIME = TIME - DEPTIME(K) 


C....CONDUCTIVE HEAT FLUX USING 2ND ORDER CURVE FIT TO RESULTS OBTAINED FROM GLENN 
WHITE'S 
C.....CODE EVAP.F 


QD = 1.4*TEMP(K)**2 + 170.0*TEMP(K) - 21300 
QI = KS*(TSO - TCPYDELTA 


ARGI = DIST/RADO 
ARG2 = (ALPHAS*LTIME)**.5/RADO 
IF (LTIME LE. LTAU(K.TEMP)) THEN 
ARG3 = 0.0 
ELSE 
ARG3 = ((ALPHAS * (LTIME - LTAU(K,TEMP)))**0.5) / RADO 
END IF 
CALL QAGI (F, BOUND, INF, EPSABS, EPSREL, RESULT, ABSERR, NEVAL, IER, LIMIT, LENW, LAST, 
TWORK, 
WORK) 
U =-0.9 * (QD + QI) * RADO * RESULT / KS 


RETURN 
END 


SUBROUTINE FAR(TIME, K, DIST, DEPTIME, U, TEMP) 
C..... THIS SUBROUTINE CALCULATES THE TEMPERATURE DEPRESSION AT A POINT DUE TO A DROPLET 
C...... LOCATED 
C.....GREATER THAN 10 DROPLET RADIT AWAY. 
C..... THE SOLUTION IS THE CLOSED FORM SOLUTION GIVEN BY CARSLAW AND JAEGER FOR A POINT 


C.....SINK. 
COMMON TDROP, RADO, TSO 
COMMON /CONST/ ALPHAS, CS, DELTA, KS, LAMDA, PL RHOS, RHOW, TCP, V0 


REAL TIME, LTIME, DIST, U, DEPTIME(1000), TEMP(1000) 

INTEGER K, TDROP 

REAL V0, RHOW, RHOS, LAMDA, CS, ALPHAS, Q, QD, QL DELTA, KS 
REAL TCP, TSO 

REAL A, B, PL, RADO 


REAL LTAU 
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LTIME = TIME - DEPTIME(K) 


IF (LTIME LT. 0.6*LTAU(K,TEMP)) THEN 
U=00 

ELSE 
QI = KS * (TSO - TCPYDELTA 
QD = 1.4*TEMP(K)**2 + 170.0*TEMP{K) - 21300.0 
Q= (QD+Q) * PI * LTAU(K, TEMP) * (RADO)*#2 
A = Q(4*RHOS*CS*(PI * ALPHAS * (LTIME -.6*LTAU(K,TEMP)))**1.5) 
B = (DIST**2)/(4 * ALPHAS * (LTIME - .6*LTAU(K,TEMP))) 
U=-A * EXP(-B) 

END IF 

RETURN 

END 


BLOCK DATA 
REAL ALPHAS, CS, DELTA, KS, LAMDA, PL RHOS, RHOW, TCP, V0 
COMMON /CONST/ ALPHAS, CS, DELTA, KS, LAMDA, PI, RHOS, RHOW, TCP, VO 
DATA ALPHAS, CS, DELTA, KS, LAMDA, PL RHOS, RHOW, TCP,. V0 /5.79E-07, 888.9, 0.0254, 
+1.297, 2225000., 3.141592654, 2520., 998.2, 35., 9e-09/ 


END 
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latter parameter is found to have a negligible effect on the droplet's evaporation. The 
droplet shape model and a modified radiation heat term are incorporated into a previously- 
developed computer model to predict the evaporation of a single droplet on a semi-infinite 
solid subjected to radiant heat input. The code predicts transient temperature profiles 
that agree well with experiment. A simplified, closed-form solution for the prediction of 
surface temperatures in the vicinity of an evaporating droplet is used to fit the data 
produced by the single droplet model. This closed-form solution facilitates calculations 
used in a model for the cooling of a surface by the evaporation of a sparse spray of water. 
The previously collected data base for sparse spray cooling using degassed water is expanded 
with a new set of experiments using water containing dissolved gases. 
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